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Supernova theory, numerical and analytic, has made remarkable progress in the past decade. 
This progress was made possible by more sophisticated simulation tools, especially for neutrino 
transport, improved microphysics, and deeper insights into the role of hydrodynamic instabilities. 
Violent, large-scale nonradial mass motions are generic in supernova cores. The neutrino-heating 
mechanism, aided by nonradial flows, drives explosions, albeit low-energy ones, of ONeMg-core and 
some Fe-core progenitors. The characteristics of the neutrino emission from new-born neutron stars 
were revised, new features of the gravitational-wave signals were discovered, our notion of supernova 
nucleosynthesis was shattered, and our understanding of pulsar kicks and explosion asymmetries was 
significantly improved. But simulations also suggest that neutrino-powered explosions might not 
explain the most energetic supernovae and hypernovae, which seem to demand magnetorotational 
driving. Now that modeling is being advanced from two to three dimensions, more realism, new 
perspectives, and hopefully answers to long-standing questions are coming into reach. 
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I. INTRODUCTION: ROOTS AND QUESTIONS 

When, why, and how can the catastrophic infall of the core of a massive star be reversed to trigger the powerful 
ejection of the stellar mantle and envelope in a supernova (SN) explosion? This fundamental problem of stellar 
astrophysics has been a matter of intense research since the crucial role of SNe for the synthesis of heavy elements and 
for the dissemination of the nuclear burning products of stars had been recognized by Burbidge et al. [lj. The latter 
authors also noticed that nuclear statistical equilibrium in the hot, dense core of evolved stars (at T > 1 x 10 9 K) 
favors iron dissociation to alpha particles, and they concluded that the huge demand of energy (about 1.7MeV per 
nucleon or 1.7 x 10 18 erg per gram) must be supplied by gravitational binding energy, leading to a contraction of the 
stellar core and ultimately to a dynamical implosion on a timescale of less than a second, f co ii ~ 0.21/^/psS, when 
the average density p 8 = p/(10 8 g/cm 3 ) exceeds unity. This groundbreaking insight is in line with Baade & Zwicky’s 
earlier idea that SNe could represent the transition of ordinary stars to neutron stars (NSs) [3|. 

Already in 1960 Hoyle & Fowler [J| proposed the two basic scenarios of stellar death: thermonuclear runaway at 
degenerate conditions (which, as we know now, drives the destruction of white dwarf stars in Type la SNe) and the 
implosion of stellar cores (associated with what is called core-collapse supernovae (CCSNe) of Types II, Ib/c, and 
hypernovae 1 ). They hypothesized (following 0 ) that the gravitational compression of the core raises the temperature 
such that thermonuclear fuel could be ignited to release the energy for triggering the ejection of the outer parts of the 
star. They also mentioned simulations by Colgate & Johnson (4), : 5|, in which the “bounce” of a forming NS launched 
a spherical shock wave that reversed the infall of the overlying stellar shells to make them gravitationally unbound. 
Colgate & White [6] realized that gravitational binding energy of order E\> ~ GM^ S /R ns > 10 53 erg, which is released 
when the core of a star collapses to a NS, is converted to neutrino emission and provides a huge energy reservoir for 
powering the SN blast wave. They argued correctly that in stellar layers pulled inward at supersonic speed along with 
the imploding core, thermonuclear combustion is unable to initiate an outward acceleration. Instead they proposed 
that a fraction of the intense neutrino flux may get absorbed in the mantle of the star to cause the explosion. 

More than four decades of theoretical and numerical modeling work, spearheaded by early pioneers of the field like 
Dave Arnett, Jim Wilson, Hans Betlre, Gerry Brown, Steve Bruenn, Wolfgang Hillebrandt, Jim Lattimer, and David 
Schramm, have helped to sharpen our picture of the diverse physical ingredients and processes that play a role in 
the core of dying stars, among them magnetohydrodynamic (MHD) effects, fluid instabilities and turbulent flows, the 
finite-temperature equation of state (EoS) of NS matter, neutrino transport and neutrino-matter interactions, and 


1 Observationally, SNe II exhibit strong H-Balmer lines in their early spectra, whereas SNe I show no H-lines. In SNe la there are Si-lines, 
in SNe lb no Si- but He-lines, and in SNe Ic none of these, indicating explosions of stars that had lost their hydrogen envelope or both 
the outer hydrogen and helium shells before collapse. More sub-classes have been introduced, some of them motivated only by recent 
discoveries: SNe II-P and II-L are discriminated by a plateau phase or linear decay of their lightcurves after the peak, lib events have 
only thin H-shells left, and spectra of Ha and Iln cases possess signatures of a dense circumstellar medium. 
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general relativistic gravity. While the bounce-shock mechanism is not supported by any modern simulation with state- 
of-the-art treatment of the physics, the “delayed neutrino-heating mechanism” as discussed by Bethe and Wilson [?} 
and aided by violent, nonradial mass motions in the collapsing stellar core (MU , has advanced to the widely favored 
scenario for powering the majority of SNe. 

The momentum behind the quest for solving the puzzle of the SN mechanism originates from important questions 
at the interface of astrophysics and nuclear, particle, and gravitational physics, for example: 

• What is the link between the properties of SNe and their progenitor stars? 

• Which stars collapse to black holes (BHs) instead of NSs, which fraction of stellar collapses do not yield explo¬ 
sions? 

• What are the birth properties of the compact remnants, i.e. their masses, spins, magnetic fields, and recoil 
velocities? 

• How can the high velocities of young pulsars be explained? Is any exotic physics necessary? 

• What characteristics does the neutrino burst from a SN have and what does it tell us about neutrino properties 
and the extreme conditions in the newly formed NS? 

• What is the gravitational-wave signature of a stellar collapse event and which information can we extract about 
the dynamical processes in the SN core? 

• What is the nucleosynthetic role of massive star explosions in the chemogalactic history? 

• Are SNe the long-sought sources of r-process elements, in particular also of the lanthanides, the third abundance 
peak, and actinides? 

• What is the population-integrated energetic footprint left by SN explosions in the dynamical evolution of galax¬ 
ies? 

In the following sections we will review the known types of stellar collapse events (Sect. |TT]), the ingredients and 
current status of numerical modeling (Sect. Mil) , the mechanisms by which massive stars might explode (Sect. IIVI) . 
and the signatures of the explosion mechanism that might serve for observational diagnostics 1 Sects. IVl and fVTl) . We 
will provide an update of recent developments as follow-up and supplement of previous reports that have approached 
the topic from different perspectives |12|-ll9i|. 

II. ROUTES TO STELLAR CORE COLLAPSE 

Massive stars possess finite lifetimes of millions to tens of millions of years, which are mainly determined by the 
period the star spends on the main sequence (MS) during central hydrostatic hydrogen burning. The evolution time 
of stars scales approximately like t evo i ~ 7.3 x 10 9 yr (M*/Mq)/(L*/Lq) with the stellar mass M* and luminosity 
L*/Lq ss (M*/Mq ) 3 ' 5 (where Mq = 1.989 x 10 33 g and Lq = 3.85 x 10 33 erg/s are solar mass and luminosity). When 
hydrogen in the stellar core gets exhausted and the star leaves the MS, its evolution speeds up considerably because 
the efficiency of energy production in the higher stages of nuclear burning decreases and concurrently energy losses 
through neutrino-antineutrino pairs rise dramatically. This is the case in particular when the central temperature of 
the star climbs to T c ~ 10 9 K, at which time e + e _ pairs become abundant and the energy drain in vv pairs accelerates 
with T 9 . At this time neutrino losses exceed the radiation losses of the star and the evolution of the helium core 
decouples from that of the stellar envelope. 

The energy drain happens at the expense of gravitational binding, leading to continuous contraction of the stellar 
core, which is slowed down only temporarily by the periods of nuclear burning. As long as nondegenerate particles 
dominate the pressure of the stellar plasma, hydrostatic equilibrium requires that the central temperature, T c , and 
central density, p c , roughly follow the proportionality 

rji 3 

— oc M? ~ const . (1) 

Pc 

According to this relation more massive stars with bigger He-cores (larger M c ) are hotter (Fig. [[]). For sufficiently 
high central temperature, nuclear fuel can ignite in the next burning stage, building up heavier and more stable 


3 



log p c [g/cm 3 ] 

FIG. 1: Stellar death regions with schematic stellar evolution tracks in the plane of central density (p c ) and central temperature 
(T c ). Colored death regions are labeled by the instability process causing the collapse of the stellar core, and the blue tracks are 
labeled by the corresponding rough birth-mass range of objects reaching the different stages of central burning (indicated by 
red dashed lines). Yellow diagonal lines mark the beginning of degeneracy (short-dashed) and strong degeneracy (long-dashed) 
of the electron plasma. Note that realistic stellar tracks exhibit wiggles and loops when the ignition of the next burning stage 
is reached and the stellar core adjusts to the new energy source (see Ref. (2pil.j 


elements in their inner core. If, however, the stellar interior enters the regime of electron degeneracy before 2 (yellow, 
short-dashed line in Fig. [[]) it ends as a white dwarf, being stabilized by lepton degeneracy pressure and cooling at 
essentially fixed density. 

Stars beyond certain birth-mass limits can reach the “death zones” in the upper and right parts of Fig. |T) where 
the stellar core becomes gravitationally unstable. Contraction, and in the case of a runaway process finally collapse, 
sets in when the effective adiabatic index drops below the critical value of 4/3 for mechanical stability (the actual 
value is slightly decreased by rotation and increased by general relativistic gravity). 

Three different processes can initiate the implosion of stellar cores in three areas of the p c -T c -plane indicated by 
different colors in Fig. [I] playing a role in different kinds of CC events. 


2 Fermions approach the degeneracy when their Fermi energy begins to exceed the thermal energy k^T. i.e. at Tg ~ 4 for nonrelativistic 
electrons and at Tiq ~ Pg / 3 f° r relativistic ones with T x = T /( 1 0 X K) and p y = ft/ { 1 (Y y gem -3 ). 
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FIG. 2: Core-density profiles of different SN progenitors at the onset of gravitational collapse. The black line corresponds to 
the ONeMg core of an 8.8 Mq star [21], the other three are SN progenitors with iron cores: an 8.1 Mq ultra metal-poor (10 -4 
solar metallicity) star (A. Heger, private communication) and 11.2 Mq fz| and 15 M© [23] solar-metallicity stars. The steps 
and kinks in the curves correspond to composition-shell interfaces (Fe/Si and O/C for the 11.2 and 15 Mq models and inner 
and outer boundaries of a C-O-Ne-layer for the 8.1 M© case). 


A. Electron-Capture Supernovae 

The lowest-mass progenitors of CCSNe develop oxygen-neon-magnesium (ONeMg) cores through C-burning fUl [H, 
m but reach e-degeneracy before hydrostatic Ne-burning can be ignited. Due to low reaction thresholds of Ne and 
Mg, the increasing electron Fermi energy enables e-captures (right upper corner of Fig. |TJ), triggering gravitational 
collapse and resulting in an electron-capture SN (ECSN). Solar-metallicity stars 3 with 9-9.25 Mq are estimated to 
have that destiny [29(, but the mass window is expected to shift and widen for lower metallicities 30] and in binary 
systems with mass loss or transfer [3l| so that ECSNe could contribute even 20-30% of all SNe [32 . l33j. 

Because of the extremely steep density decline in a thin C-O-shell (~0.1 Mq between about 3 x 10 4 g/cm 3 and 
4 x 10 8 g/cm 3 ) at the edge of the O-Ne core (Fig. [2]), these stars have special explosion properties (Sect. HVCl) . They 
eject little carbon and oxygen and very little nickel, their SNe will therefore be relatively faint. The Crab remnant of 
SN 1054 is thought to be the relic of such an explosion [HI HI], and an increasing number of dim events like SN 1997D, 
1999br, 2005cs, 2008S and other recently observed transient sources are discovered as possible candidates. 
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The metallicity Z is the total mass fraction of chemical elements heavier than helium in the matter the star was formed of. The solar 
value has been determined to be 0 . 016 . 
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Main Sequence Mass (M 0 ) 




FIG. 3: Kinetic energies and ejected nickel masses for stellar explosions versus initial (zero-age main sequence; ZAMS) masses 
from different authors. While explosion properties are deduced from comparing observations with lightcurve and spectra 
calculations based on (spherically symmetric) models, the ZAMS masses are estimated by linking ejecta masses to initial 
masses through stellar evolution models with mass-loss assumptions (upper panels [HI and lower left panel [Iff]; reproduced 
by permission of the authors and AAS and ©ESO, respectively) or by inferring ZAMS masses or upper limits from computed 
stellar evolution histories that account for the properties of discovered SN progenitors or their stellar environments (i.e., coeval 
star clusters, host galaxies [2<|,|22}; lower right panel, kindly provided by John Eldridge and Stephen Smartt). Due to theoretical 
uncertainties this leads to largely different mass determinations for some cases (SN 1987A, SN 1999br, SN 1999em, SN 2004et, 
SN 2005cs). For masses >25 Mq Tanaka et al. [H] discriminated a very energetic and bright HN branch from a low-energy, 
faint SN branch. The objects with little nickel production of the latter branch, however, have also been interpreted as weak 
explosions (possibly ECSNe) near the lower mass limit for SN progenitors. 


B. Iron-Core Supernovae 

Massive stars that ignite hydrostatic Ne-burning form an iron core. The latter becomes gravitational unstable when 
the nuclear statistical equilibrium (NSE) at temperatures around 10 10 K (k^T ~ 1 MeV) favors the dissociation of 
iron-group nuclei to a-particles and a growing number of free nucleons (upper region in Fig. [[]). With the onset of 
contraction and increasing density and electron chemical potential, also e-captures on nuclei (and some free protons) 
speed up and accelerate the implosion. The dynamical collapse is abruptly stopped only when nuclear densities 
(p > 2.7 x 10 14 g/cm 3 ) are reached, and the phase transition to homogeneous nuclear matter leads to a sudden 
increase of the effective adiabatic index due to repulsive short-range forces between nucleons. 

When the overshooting inner core rebounds and crashes supersonically into the subsequently infalling layers, sound 
waves steepen into a shock front that ultimately leads to the disruption of the star in the SN explosion. However, 
different from ONeMg cores, the much flatter density profile in and around Fe cores (Fig. [2]) leads to long-lasting, 
high mass accretion rates and large ram pressure of the infalling shells. This impedes the outward propagation of the 
shock and makes Fe-core progenitors harder to blow up than stars with ONeMg core. Although more massive stars 
exhibit a gross tendency to bigger He-cores and shallower density decline, the variation with stellar birth mass is not 
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necessarily monotonic [22J. The mechanism(s) by which ECSNe and Fe-core SNe succeed to explode will be discussed 
in Sect. UVl 

Stellar cores of pre-SN stars are expected to rotate relatively slowly, i.e. with average pre-collapse spin periods 
of tens of seconds or more. This is a consequence of angular momentum loss associated with mass loss phases (in 
particular when the star becomes a red giant), because magnetic torques from fields generated by differential rotation 
in the star couple core and envelope and thus transport angular momentum efficiently out of the core [36| . Stellar 
rotation is therefore not expected to play a crucial role for the explosion mechanism of normal CCSNe (Sect. IIVDI) . 


C. Gamma-Ray Burst Supernovae 

Rapid stellar rotation, however, is thought to be crucial in the case of gamma-ray burst (GRB) SNe and hypernovae 
(HNe) (for a review, se e (1511 1. The latter originally obtained their name because of exceptional brightness and thus 
high nickel production [37| but are now considered as stellar explosions with unusually high ejecta velocities (i.e., 
very broad spectral lines) and thus large kinetic energies [H( (Fig- 0)- HNe are found to be associated with long- 
duration (tGRB > 2 s) GRBs, either observed spectroscopically (e.g., SN 1998bw with GRB 980425, SN 2003dh with 
GRB 030329, SN 20031w with GRB 031203, SN 2006aj with GRB 060218, SN 2010bh with GRB 100316D) or as late 
lightcurve humps superimposed on the power-law decline of the afterglow that follows the GRB. 

GRBs with their extremely luminous high-energy radiation are understood as ultrarelativistic, collimated outflows 
(“jets”). Strong global asymmetry is also suggested by line profiles (in particular double-peaked oxygen emission 
lines) seen in many HNe. Such events are interpreted as signatures of BH-forming stellar collapses (“collapsars” ®), 
in which matter around a rapidly spinning BH is able to set free energy in neutrinos, electromagnetic poynting 
flux, and mass outflow with an efficiency of up to roughly 40% of the rest-mass energy of accreted material, E acc < 
0.4Mc 2 ~ lO 54 (Af/M 0 /s) erg/s. Alternatively, a nearly critically rotating NS, t ro t ~ 1ms, with ultrastrong dynamo¬ 
generated magnetic field, (B) > 10 15 G, (“millisecond magnetar”) is discussed as possible central engine of GRBs and 
HNe. The jet and stellar explosion could either be powered by rotational energy of the magnetar or by gravitational 
and rotational energy of the accretion flow and BH. Both can be tapped by magnetic fields through MHD effects 
(Sect. HVDl) and by neutrinos radiated from matter heated by magnetically generated viscous dissipation j4Cj. The 
existence of expected strong disk “winds” with the observed large Ni production [l5| , however, seems to be challenged 
by MHD simulations fill ]. 

The progenitors of collapsars and GRB-HNe are thought to possess a massive core, which forms a BH instead 
of exploding before. They must be compact stars without an extended hydrogen envelope in order to allow jets 
to emerge ultrarelativistically, i.e., the crossing time of the jet must be shorter than the on-time of the central 
engine: i?*/c < f e ngine- Moreover, the collapsing stellar core must contain a high specific angular momentum, 
j > GMbh/c > lO 16 M BH /(3M 0 )cm 2 /s, to either form a magnetar with the necessary huge reservoir of rotational 
energy or to allow for a thick, massive accretion disk that remains long enough around the newly formed BH to 
efficiently release energy. 

Such requirements favor rapidly rotating Wolf-Rayet stars as progenitors, but special initial conditions (a high birth 
spin) and evolution paths avoiding combined mass and angular momentum loss, or alternatively binary scenarios, are 
necessary [Z^, EH- In the present-day universe HNe and GRB-SNe are rare (with a GRB/SN ratio of ~1/1000, less 
than 1% of all SN Ib/c produce GRBs), but BH-forming CC events and GRBs could be very common in the early 
(metallicity less than ~1/10 solar) universe. This theoretical expectation is compatible with the fact that GRB-SNe 
are preferentially (but not exclusively) observed in low-metallicity environments. 


D. Pair-Instability Supernovae 

Stars above about 100 M 0 are very hot and encounter the pair instability (in the upper left corner of Fig. [[]) after 
central carbon burning (e.g., [ 22 ], 0] and refs, therein) at T ~ 10 9 K. The gravitational instability occurs because the 
formation of e + e _ pairs from high-energy photons converts thermal energy to rest-mass energy and thus reduces the 
adiabatic index of the EoS below 4/3. 

In the stellar mass range between ^100 M 0 and ~14OM 0 and for M* > 260 M@ collapse to a BH is expected. For 
intermediate masses the ignition of the still available thermonuclear fuel during the implosion is violent enough to 
trigger the complete disruption of the star with an explosion energy up to more than 10 53 erg and the production of 
up to >50 Mq of 56 Ni f2p[ [43] . While such “thermonuclear core-collapse SNe” were originally termed “hypernovae” 
by Woosley & Weaver [4h|, they are now commonly called pair-instability SNe (PISNe) or pair-capture SNe (PCSNe). 
In the case of BH formation, in particular in the presence of rotation allowing for an accretion torus, huge amounts 
of energy are released in neutrinos, roughly (0.01-0.03)M*c 2 ~ 10 55 erg, depending on the angular momentum fl(| . 
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TABLE I: Neutrino reactions with stellar-medium particles and between neutrinos in the Garching models. N means either n 
or p, v € {is e , t'e, ty,, v T , ty}, and v x € -fry, , v T , u T }. In addition to “inelastic” nucleon recoil, thermal motions, phase-space 
blocking, high-density A-77-correlations [5f] and weak magnetism corrections [53 ]■ also quenching of the axial-vector coupling 
[53| and the reduction of the effective nucleon mass at high densities f54| are taken into account in the rates marked with a 
dagger (*). A prime indicates that the neutrino can exchange energy with the scattering target (non-conservative or “inelastic” 
scattering) 


Process 

References 

Beta-Processes 

t'e + n ^ e~ + p 

[53J + 

u e +p ^ e + +n 

[51] f 

v e + (v4, Z ) ^ e + (A, Z + 1) 

[55] 

Scattering Reactions 
v + (v4, Z ) ^ v* + (A, Z ) 

[56] (ion-ion correlations) 

v + N ^ u' + N 

[57] (inelastic contribution) 

[51]+ 

^ + e ± ^ v' + e* 

[58] 

(“Thermal”) Pair Production 
v + v ^ e~ + e + 

[59, 60] 

Nucleon-Nucleon Bremsstrahlung 
v + v + N + N ^ N + N 

m 

Reactions between Neutrinos 

m 

Vx+{v e ,V e } ^ u' x + {v’ e ,K} 

[62] 


While some recently discovered ultra-bright SNe and transients (for example SN 2002ic, 2005gj, 2005ap, 2006gy, 
2007bi, 2008es, 2010gx) have been discussed as PISN candidates (e.g., [43, |48}), other explanations for the extreme 
luminosity than excessive Ni yields have been proposed, e.g., interaction of explosion ejecta with a dense circumstellar 
medium [49j] or additional energy release by magnetar spin-down f5(il ]. The expected rate of PISNe is small, maybe 
one of 100-1000 normal stellar core collapses, and presumably mostly associated with metal-poor host galaxies. In 
the Milky Way possibly two dozen very massive hypergiants like the evolved luminous blue variable star rj Carinae 
might end their lives in such events. 

Although the nature of the stellar death events associated with ultra-bright transients and in particular the energy 
source of their extraordinary luminosity are not at all clear and will remain a topic of intense research and debate in 
the coming years, space limitations demand to constrain the rest of this article mostly on the physics and processes 
that are relevant for the far majority of ordinary CCSNe. 


III. NUMERICAL MODELING AND PHYSICS INGREDIENTS 


While over two decades that followed the poineering work by Colgate & White [y , Arnett 
modeling was constrained to spherically symmetric (ID) simulations with few exceptions only 
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and Wilson [g4|, SN 
6a l, the situation has 


radically changed in the post-SN 1987A era. Detailed observations of this nearest SN in the era of modern astronomy 
revealed that large-scale mixing processes had transported radioactive nuclei with velocities up to ~4000km/s from 
the deep core far into the hydrogen envelope of the exploding star, suggesting that spherical symmetry was broken 
already during the very first moments of the blast [69j,[70|. Moreover, two-dimensional (2D) simulations in the early 
1990s demonstrated that violent convective overturn takes place in the neutrino-heating layer between the gain radius 
and stalled accretion shock mm. This raised hopes that buoyant energy transport to the shock could crucially 
support the delayed neutrino-heating mechanism and finally ensure robust explosions after ID models had turned 
out to be successful only with special assumptions that could not withstand closer and more detailed analysis. For 
example, neutron-finger instability inside the nascent NS was proposed to enhance the neutrino luminosities and thus 
neutrino heating 0 but is disfavored because lepton equilibration between fingers and surroundings was shown to 
proceed faster than thermal equilibration i& 

In the following sections more recent developments and the present status of numerical approaches will be briefly 
summarized. While three-dimensional (3D), general relativistic (magneto-)hydrodynamic simulations including mi- 









crophysical EoS and sophisticated, energy-dependent neutrino transport are the ultimate, brave objective, only first 
steps have so far been achieved, approaching the goal from different directions. Mastering this grand computational 
challenge will require highly parallelized codes with excellent scaling capability on tens of thousands of processor cores 
to achieve sustained performance on the hundreds of terafiop/s to petaflop/s level. Still, one 3D model calculation 
will take a wall-clock time of several weeks to months. 


A. Hydrodynamics and Gravity 


To date fully self-consistent modeling of stellar collapse and explosion in 3D has been achieved only by Fryer 
and collaborators [73i475l |. yet only by sacrificing many aspects which are important for quantitatively reliable and 
conclusive results concerning the SN mechanism. In particular, Newtonian gravity and a grey, flux-limited neutrino 
diffusion (FLD) scheme JJ were applied in combination with a smoothed particle hydrodynamics (SPH) method, 
which permits economical calculations in 3D with relatively low resolution but is noisy and diffusive. Good resolution 
and an accurate representation of the hydrodynamical quantities, however, are essential to treat the growth of fluid 
instabilities from initial seeds in the SN core [zi[z3, and general relativity as well as a multi-group description of 
neutrino transport including velocity-dependent observer corrections were found to cause important differences in 
ID [78;, 79j and 2D simulations ;n[). 8~fl . 

Other groups, using mesh-based discretization schemes for solving the hydrodynamics, have so far studied only 
more constrained problems in 3D than Mzl, mostly also making even more radical approximations of the relevant 
(micro)physics. For example, in refs. 82, 83jjthe development of a nonradial hydrodynamic instability of the accretion 
shock in a collapsing stellar core, the so-called standing accretion shock instability (SASI; 0),' was investigated 
for a steady-state flow thr oug h outer and inner grid boundaries with an ideal-gas EoS and parametrized neutrino¬ 
cooling terms. In refs. j85l I86| a similar accretion setup was studied with a microphysical EoS and additional simple 
neutrino-heating terms for prescribed luminosities and spectra (neutrino “lightbulb” approximation, NLA), which 
enabled neutrino-driven convection. Using the NLA, refs. 0,113 investigated the onset of an explosion in “realistic” 
collapsing stellar cores systematically by varying the driving neutrino luminosity to explore the dependence on the 
dimension (ID, 2D or 3D) of the simulation. Gravitational-wave (GW) signals from infall, core bounce, and early 
postbounce (~400 ms) phases were computed with 3D general relativity (GR) for NS and BH formation [1301 an d 
with 3D Newtonian^hydrodynamics and an effective general relativistic potential (developed as an approximation of 
GR gravity in 0, [92]) for NS formation [§3}, making various crude simplifications of the neutrino effects and partly 
even of the EoS of the stellar plasma. 

A grey description of the neutrino transport outside of an excised high-density core of the proto-neutron star (PNS) 
according to ref. [94] was applied for exploring NS kicks, neutrino emission asymmetries, and GW signal characteristics 
by long-time 3D simulations of SN explosions in refs. 0,[96|. Very first results of Newtonian 3D calculations with more 
detailed multi-group (MG) transport treatments have already been put out, using “ray-by-ray” (RbR) MGFLD 97] 
or an implementation of the “isotropic diffusion source approximation” (IDSA, see Sect. IIILBl) for v e and v e with 
or without [99] a RbR approach, coupled to a trapping treatment for heavy-lepton neutrinos. 


B. Neutrino Transport 


Over the past decade sophisticated multi-energy group solvers for three-flavor neutrino transport including energy- 
bin coupling and velocity-dependent terms (corrections due to the motion of the stellar plasma) have been developed 
and applied to all stages of stellar core collapse and the transition to explosion in ID calculations. On the one hand this 
was achieved by direc t int egration of the Boltzmann transport equation (BTE) with a discrete-ordinate (Sn) method in 
GR simulations [TfirS L fioi| . on the other hand by integrating the set of two-moment equations of the BTE for neutrino 
number, energy, and momentum using a variable Eddington-factor closure obtained from convergent it erat ion with 
a model (i.e. simpli fied ) Boltzmann equation. The latter approach was developed for Newtonian 0, Il02| ] as well 
as GR simulations |l03j . It was also generalized for multi-dimensional applications by adopting a “ray-by-ray plus” 
(RbR+) approximation 0,0, in which spherical transport problems are solved on each angular bin of a 2D or 3D 
polar coordinate grid. This approximation implies that the neutrino intensity is assumed to be axially symmetric 
around the radial direction and the neutrino flux is considered to be purely radial. The “plus” suffix signals, however, 
that neutrino pressure gradients and the lateral advection of neutrinos with fluid flows are taken into account in the 
optically thick regime to prevent artificial hydrodynamic instabilities 0|. 

All the published ID and 2D SN models of the Garching group, e.g. in refs. SEoira , include the full, state- 
of-the-art set of neutrino interactions listed in Table Q] Recently, ID results based on a similarly refined treatment of 
the neutrino processes have been put out by the Oak Ridge group [79]. 
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Truely multi-dimensional, energy-dependent transport schemes for radiation-hydrodynamics with neutrinos have so 
far been used extensively in 2D Newtonian simulations only by the Arizona-Hebrew-Princeton collaboration, which 
applie d a MGF LD method (e.g., in refs. and an Sn solver for a multi-angle (MA) treatment treat¬ 

ment uni, [Hi, however without energy-bin coupling and without properly accounting for effects associated with 
fluid motions. This are severe shortcomings |79l |80||. which are avoided in more elaborate 2D Newtonian implemen¬ 
tations of MGF LD |l!3| and of a two-moment closure scheme for the coupled set of neutrino energy and momentum 
equations |ll4j . An alternative approach is the “isotropic diffusion source approximation” (IDSA; [l!5l |). in which 
the neutrino distribution function is decomposed into trapped and streaming particle components, whose separate 
evolution equations are coupled by a diffusion source term. This method was simplified to a RbR version for v e and 
v e and only a subset of neutrino processes in 2D [116} and 3D [98} SN simulations. A more detailed comparison and 
critical assessment of presently employed transport treatments can be found in ref. (79|j. 

Pointing into the future, routes towards 3D time-dependent neutrino transport in radiation-hydrodynamics cal¬ 
culations have been outlined in the form of a rigorous solution of the 6+1 dimensional (three spatial dimensions, 
energy and two direction angles for the radiation momentum, plus time) BTE by an Sn discretization scheme EH, 
by spectral methods LUS}, and in GR by a truncated moment formalism uM , but observer corrections due to fluid 
motion, relativistic effects, nonlinear, energy-coupling interaction kernels, and high parallization efRcieny are major 
challenges. 

Direct comparisons of multi-dimensional SN calculations by different groups with different codes and approximations 
have not been carried out so far, in contrast to the ID case fioilml. and will be a formidable task for the coming 
years. However, FLD was shown to underestimate angula r var iations of the radiated neutrinos and to sphericize the 
radiation field compared to a multi-angle (Sn) treatment m, although fundamental changes of the hydrodynamic 
evolution were not observed despite higher neutrino-heating rates with the Sn code. On the contrary, the RbR 
approximation generically sharpens angular variations since all fluxes are radial. Local emission maxima (“hot spots”) 
in the neutrinosplieric region therefore send radiation only in radial direction. Nevertheless, because nonspherical 
accretion flows in the SN core exhibit unsteady behavior in space and time (see, e.g., [9(|) “variational averaging” 
can be expected to diminish any dynamical consequences of local emission peaks [80j] ■ 

Despite undeniable weaknesses, the complexity and computational intensity of neutrino-hydrodynamics in full 
generality will make the use of simplifications unavoidable still for some time. 


C. Equation of State and Composition of Stellar Plasma 

The nuclear and subnuclear EoS is an extremely important ingredient for SN modeling. Unfortunately, our knowl¬ 
edge of, in particular, the supranuclear regime is incomplete, although information from nuclear theory and experi¬ 
ment [l23] as well as astrophysical observations, for example by the recently discovered 1.97 Mq binary millisecond 
pulsar PSR J1614-2230 Il24j . is rapidly growing and is beginning to set serious constraints on the possible existence 
of bigger mass regions with exotic phases in NS interiors [l25l ll26| . 

The two EoS s for hot NS mat ter w idely used for stellar core collapse in the past decade are those of Lattimer 
& Swesty and Shen et al. |l28} . They include nuc leons and nuclei, electrons and positrons, and photons. 
The former is based on a compressible liquid-drop model [l29t | with a Skyrme force for nucleon interactions. The 
transition to homogeneous nuclear matter is established by a Maxwell construction. Most of the SN simulations by 
the Garching group were performed with a version (LSI8O-E0S) with an incompressibility modulus of bulk nuclear 
matter of K = 180 MeV and a value of 29.3 MeV for the symmetry energy parameter. The Shen et al. EoS (STOS- 
EoS) employs a relativistic mean field model with parameter settings that reproduce characteristic properties of heavy 
nuclei. It is extended with the Thomas-Fermi spherical-cell approximation to describe homogeneous matter as well 
as inhomogeneous conditions. Its bulk incompressibility and symmetry energy have values of 281 MeV and 36.9 MeV, 
respectively. 

These EoSs describe the nuclear composition as a mix of free nucleons, alpha particles, and a representative heavy 
nucleus, whose mass and charge numbers (A, Z) depend on density, temperature, and neutronization of the matter. 
Although largely different (A, Z) are returned by both EoSs during the infall stage and affect, for example neutrino 
trapping through coherent neutrino-nuclei scatterings, ID simulations yield basically the same behavior. Quantitative 
differences occur only on a modest level of 5-25% in quantities characterizing collapse, bounce, and early postbounce 
evolution, e.g. in the central lepton fraction at neutrino trapping, position of shock fo rmat ion, peak luminosity of the 
v e burst, and maximum radius to which the shock expands before it retreats again Hamas This outcome is even 
more astonishing in view of the appreciably different adiabatic index T = (dhiP/dhip) s (P, p 1 and s are pressure, 
density, and entropy per nucleon, respectively) for both EoSs around nuclear density (Tls ~ 2.2, Tstos ~ 2.9) and 
the correspondingly different maximum compression and rebound behavior at bounce. Once again Mazurek's Law 
applies, according to which any change of the microphysics is moderated in its effects on collapsing stellar cores by a 








10 



2000 1000 a 10Q0 2000 2B0B 1000 0 1000 2000 

r [km! r [km! 





FIG. 4: Two-dimensional SN simulations Il2ll1 of an 11.2 Mg star 22] for three different nuclear EoSs. The upper panels 
show cross-sectional entropy distributions at 412 ms after b ounc e for the LSI 8 O-E 0 S (left), at 586 ms p.b. for the STOS-EoS 
(middle), and at 500 ms p.b. for the Hillebrandt & Wolff EoS [l22l ]. The last is the stiffest EoS of the set. It leads to the slowest 
contraction of the PNS (bottom left) and because of weaker neutrino heating and less vigorous hydrodynamic mass motions 
does not yield an explosion within the simulated time as visible in the evolution of the average shock radius (bottom right). 


strong feedback between the EoS, weak interactions, neutrino transport, and hydrodynamics 113,1 

In order to achieve a more elaborate treatment of the nuclear composition in the shock-heated region below neu- 
trinospheric densities after bounce and to connect smoothly to the chemical abundances of the progenitor star, the 
Garching simulations employ at p < 10 11 g/cm 3 a Boltzmann-gas NSE description with typically two dozen nuclear 
species, and in the non-NSE regime at T < 5 x 10 9 K a nuclear “flashing” treatment [91] or, alternatively available 
now, a small reaction network for nuclear burning. 

With a maximum gravitational mass of 1.83 Mq for cold NSs in weak equilibrium, the LSI 8 O-E 0 S is not compatible 
with PSR J1614-2230. Moreover, an incompressibility of K = 18 0 Me V seems in conflict with the experimentally 
favored value of K ~ 240 MeV for symmetric nuclear matter OH,(HI. While the STOS-EoS (M™ ~ 2 .22M 0 ) 
fulfills both constraints, its radius of ~15km for a 1.4M 0 NS does not match the best NS radius esti mat e from 
the currently most comprehensive evaluation of astrophysical data, R ns ~ 11-12.5km for $ = IAMq |l26| . This 


estimate overlaps with the range of ~ 10-14km deduced from theoretical considerations [1231 ]. which in turn agrees 
with a NS radius of ~12km for the LSI 8 O-E 0 S. 

The properties of cold, neutronized NSs, however, are not necessarily conclusive for the conditions in the hot 
SN-core environment. Indeed, for different versions of the LS-EoS with K = 180,220,375 MeV (the last two being 
compatible with PSR J1614-2230) the structure of hot PNSs well below the maximum mass, which is a relevant aspect 
for the early postbounce evolution of collapsing stellar cores, shows only smaller differences. Correspondingly, I D CC 


simulations with these EoS versions revealed only minor differences until hundreds of ms after bounce 131 , Il3(l ll37| . 


During the later PNS cooling phase and in particular when mass accretion brings the PNS close to the mass limit, 
differences in the stiffness and the symmetry energy of the EoS can have important consequences, e.g. for the time 
when BH formation occurs [l32l ] or for conv ective activity in the PNS and its influence on the neutrino emission [l38| . 
Moreover, 2D simulations showed EMEU that the explosion of 11.2 M 0 and 15M 0 progenitors depends sensitively 
on the radius evolution of the PNS in the first few 100 ms after bounce, i.e., the radius contraction of the PNS (in 
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FIG. 5: Left panel: Neutr ino-p owered ECSN of an 8.8 Mg star with ONeMg core [ 2 H. [28| visualized by mass-shell trajectories 
of a ID simulation (from 11051) '). The SN shock (bold, outgoing line) expands for ~50ms as accretion shock (the downstream 
velocities are negative) before it accelerates by reaching the steep density gradient at the edge of the core. Neutrino heating 
subsequently drives a baryonic “wind” off the PNS surface. Colored lines mark the inner boundaries of the Mg-rich layer in 
the O-Ne-Mg core (red; at ~0.72AIq), C-O shell (green; at ~1.23 Mq), and He-shell (blue; at ^1.38 Mq). The outermost 
dashed line indicates the gain radius, and the inner (bold) solid, dashed, and dash-dotted lines are the neutrinospheres of v B , 
v e , and i/ Xl respectively. Right panel: Neutrino lumi nos ities and mean energies from an ECSN for the infall, u e breakout-burst, 
accretion phase, and PNS cooling evolution (from 110711 ). The average energies are defined as the ratio of energy to number 
fluxes. (The left panel is reproduced with permission; copyright: ESO.) 


contrast to the final radius of the NS), because a more rapidly shrinking remnant radiates neutrinos with higher fluxes 
and energies mm, thus enhancing neutrino heating and in particular also enabling more violent hydrodynamic 
instabilities (Fig. [I]). 

A variety of new non-zero temperature EoSs for SN studies have been put out recently Il39l4l42j| . With modern 
ID SN codes being available, these EoSs have been (or will be) channelled through an “industrial” testing pipeline, 
confirming the modest influence of differenc es n ear or above nuclear-matter density on the early shoc k evo l utio n in 
ID as reported from previous studies above Il32j| . Also a refined description of the nuclear composition jl32lll39l ) did 
not manifest itself in a big impact on infall and shock formation. Future stu dies, also in 2D and 3D, will have to show 
whether so far ignored light nuclei ( 2 H, 3 H, 3 He, Li) besides 4 He |l43l - ll45l will have any relevant effects on the SN 
mechanism [1321 Il43| or the neutrino-driven wind from the cooling PNS m 


IV. EXPLOSION MECHANISMS 

In this section we will review the mechanisms by which the gravitational binding energy of the collapsing stellar 
core can be tapped to eject the outer stellar layers in a SN blast. A particular problem in understanding the onset of 
the explosion of massive stars is connected to the need of reversing implosion to explosion by transferring energy from 
the nascent NS to the overlying shells. This is different from thermonuclear SNe (SNe la) of white dwarfs (WDs), 
where the combustion (by deflagration or detonation) of carbon and oxygen to nickel and silicon in an essentially 
hydrostatic object releases sufficient energy to unbind and destroy the whole star. 

The typical energy scale of the explosion of a WD near its Chandrasekhar mass limit is set by the release of nuclear 
binding energy associated with the conversion of ~1 Mq C+O to Si+Ni (roughly 2 x 10 51 erg) minus the gravitational 
binding energy of the initial, highly degenerate WD (several 10 5 °erg). But what sets the energy scale of CCSNe? 
Why does the large majority of “normal” cases have explosion energies similar to SNe la? The answer to this question 
is connected to the initial state of the dying star, in which the gravitationally unstable Fe-core is a configuration 
resembling a massive, degenerate WD, surrounded by dense shells whose gravitational binding energy is of the same 
order of magnitude, i.e., around (1-15) xlO 50 erg. 

Any self-regulated mechanism for powering the explosion will deposit an energy in this range, possibly a few times 
this value, before the energy transfer is turned off. The neutrino-heating mechanism is such a self-regulated process, 
because the matter particles absorbing energy from neutrinos will react by expanding away from the heating region 
as soon as they have acquired an energy of the order of their binding energy. This will evacuate the heating region 
and diminish further energy input. But CC events show a large diversity concerning their kinetic energy, ranging 
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FIG. 6: Neutrino-driven explosions of Fe-core progenitors The upper left, upper right, and lower left panels display 

the time evolution of color-coded entropy profiles in the north and south pole directions for 2D simulations of an 8.1 Mg ultra 
metal-poor (10 -4 solar metallicity) star (A. Heger, private communication), and 11.2 Mq [2S| and 15 Mq [23:] solar-metallicity 
stars, respectively. The shock position is clearly visible as a sharp boundary between high-entropy (yellow, red) and low- 
entropy (blue, black) regions. Shock oscillations are associated with violent convective activity in the neutrino-heating region 
and strong, bipolar SASI sloshing motions of the whole postshock layer. The explosions develop highly aspherically in all cases. 
The lower right panel shows, for example, an extreme dipole asymmetry of the cross-sectional distribution of electron fraction 
(' Y e ; left) and entropy at 775 ms p.b. for the 15 Mq model, which explodes in a unipolar way. The NS is located at the position 
of the lowermost long tickmark on the vertical axis, far away from the geometrical center of the roundish shock contour (white 
line). 


from >10 5O erg to some 10 51 erg for SNe and up to several 10 52 erg for HNe (Fig. [3]). Explosion energies far in excess 
of the initial gravitational binding energy of the ejecta suggest a different driving mechanism than neutrino heating, 
a process without the self-regulation described above. Magnetorotational explosions fulfill this requirement, because 
the blast-wave energy is extracted from the huge reservoir of rotational energy of a rapidly spinning PNS by magnetic 
fields and can be delivered in the form of electromagnetic energy. 

In the following the status of our present understanding of both of these mechanisms will be summarized and also 
other, more controversial suggestions will be addressed. 


A. Thermonuclear Mechanism 


While ignition of thermonuclear combustion in compression-heated, free-falling shells was recognized not to be able 
to blow matter outwards Russian authors jl48l - ll51 later proposed neutrino radiation from the collapsing stellar 
core to heat the degenerate C+O shell of a low-mass progenitor star at hydrostatic conditions and a density around 
2 x 10 9 g/cm 3 and thus to ignite a thermonuclear burning front that explodes the star. The heating was considered 
to be mainly by neutrino-electron scattering. 

Although this is a cute idea, neither the stellar nor the dynamical conditions assumed for this scenario could be 
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verified by detailed progenitor and explosion models. In ONeMg-core progenitors, for example, which define the low- 
mass limit of stars that undergo core collapse to radiate large neutrino luminosities, the C+O shell is initially located 
between roughly 500 and 1000 km (at densities < 4 x 10 8 g/cm 3 ) and falls dynamically inward (with compression- 
induced burning) much before it is exposed to a high fluence of neutrinos (see Fig. [5]). If, in contrast, the O and C 
layers are farther out at r > 1000 km as in more massive Fe-core progenitors (Fig. [2j, the neutrino flux is diluted by 
th e larg e distance from the source and the electron densities (and degeneracy) there are much lower than adopted 


[148M151 1 so that neutrino-electron scattering cannot raise the temperature to the ignition threshold. 

Presently PISNe are the only stellar CC events where the explosion mechanism is known to be based on ther¬ 
monuclear energy release (Sect. Ill Dl) . but a closer examination of the possibility of neutrino-triggered burning in the 
significantly more compact low-metallicity stars might be interesting. 


B. Bounce-Shock Mechanism 

The purely hydrodynamical bounce-shock mechanism Si , in which the shock wave launched at the moment of 
core bounce (Sect. Ill Bl) causes the “pro mpt” ejection of stellar mantle and envelope, has been a matter of intense 
research in the 1980’s (for a revie w, se e [152j j. However, for more than 20 years now detailed analytical analysis of 
collapse and bounce physics (e.g., jl53j and refs, therein) as well as all modern CC simulations —despite continuous 
improvements and significant quantitative differences in details, mainly linked to important refinements of electron 
captures on heavy nuclei and neutrino-electron scattering during infall — agree in the basic outcome: The prompt 
mechanism cannot succeed in causing the explosion of any progenitor star. 

Upgrades of the microphysics turned out to disfavor prompt explosions by decreasing the size of the homologously 
and subsonically collapsing “inner core”, whose mass scales with the instantaneous Chandrasekhar mass, Mch{t) oc 
Y 2 (t), and whose edge defines the location of shock formation at bounce. (The number fraction 1) of particles, here of 
electrons (e), is defined as the number of particles (i) per nucleon.) With the presently most sophisticated treatment 
of neutrino emission by electron captures on nuclei and free protons during core infall [55| the central electron fraction 
Y e c after neutrino trapping decreases to 0.25-0.27 (corresponding to a trapped lepton fraction of 0.285-0.30), and 
the shock formation point (defined by the location where the entropy first reaches a value of 3 k B per nucleon) lies at 
an enclosed mass of only O.4-O.5M 0 mmm- Moreover, because preferred nuclear EoSs are relatively stiff, the 
rebound of the inner core is too weak to transfer a large energy to the shock. The flow discontinuity, running into 
supersonically infalling material at densities below ~10 13 g/cm 3 , quickly loses its initial energy by heating the plasma 
to entropies of several fee per nucleon and thus disintegrating heavy nuclei to free nucleons (which consumes roughly 
1.7 x 10 51 erg per 0.1 M 0 ). A short transient period of positive postshock velocities therefore lasts only 1-2 ms, after 
which the velocity in the whole postshock region becomes negative again. A negative postshock velocity defines the 
moment of shock stagnation, at which time the shock has travelled through just 0.3-0.35 M 0 of iron material and 
is still deep inside the stellar Fe-core. Since the preshock density then is still above 10 11 g/cm 3 , shock stagnation 
happens well before shock breakout and thus before the release of the prompt burst of v e . Lepton number and energy 
(~2xl0 51 erg) drain by the escape of the v e burst is therefore not causal for the shock stagnation. 

Despite negative velocities and thus accretion flow to the central NS in the downstream region of the shock, the 
latter continues to propagate outward in mass as well as radius. This motion of the shock stagnation radius is a 
response to the massive accretion of infalling matter (initially with a rate M ^ lM 0 /s; cf. Fig. [5]), which emits 
energy and lepton number in neutrinos and thus settles onto the PNS only gradually, creating the postshock pressure 
that pushes the shock position outwards. Finally, after reaching a maximum radius of typically 100-150 km, the 
accretion shock retreats again in ID models, following the contraction of the nascent NS roughly according to the 
relation 


(L v {$))*/*$*/* Lt /3 

S ^ M 2 / 3 Mns /3 ^ M 2 / 3 Mns /3 ^ M 2 / 3 Mnl 3 (k B T v ) s / 3 ’ 

which can be derived by combining equations (33, 39, 43, 44, 56, 63) of ref. |l55j and assuming R g oc R ns for the 
“gain radius” R g (see Sect. IIV Cl) as well as L v oc R 2 s T*, and (e 2 ) oc (k B T v ) 2 for neutrino (v G {u e ,u e }) luminosity 
and mean squared energy, respectively. The radius of maximum shock expansion at this stage is still well below the 
dissociation radius of iron, for which the equality GM ns m u /Rdiss = 8.8 MeV ( m u is the atomic mass unit) yields 
Rdiss ~ 200 km. This means that the matter behind the shock is fully disintegrated into neutrons and protons. 

It is during this period of ~100ms of slow shock expansion that a “gain radius” emerges, at which energy losses by 
neutrinos for r < R g change to neutrino heating for r > R g 0 ] ■ Before this moment neutrino losses are dominant in 
the whole postshock layer. The onset of neutrino-energy deposition also allows convective activity to develop behind 
the shock: Neutrino heating creates a negative entropy gradient ds/dr, which is unstable in the strong gravitational 
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field according to the Ledoux criterion, 


C L 


/ dp\ ds / _dp_\ d Y e 
ydsJ Ye ,pdr \dY e ) s p dr 


(3) 


This criterion, however, defines growth conditions for convection and Rayleigh-Taylor stru ctures exa ctly only for 
a static layer, but needs to be generalized for the infalling flow in the postshock region [76|, 1104 Il56| . Convective 
activity can take place there only when the inward advection timescale f a dv ~ r/|u r | for radial velocity v r is sufficiently 
longer than the convective growth timescale of perturbations (inverse Brunt-Vaisala or buoyancy frequency), t conv ~ 
(<7grav | Cl |/ p) 1 ^ 2 , or than the buoyancy acceleration timescale, fbuoy ~ (<7 g ra v (5p/p) /r) 1 ^ 2 , of blobs with density 
contrast Sp/p in a local gravitational field with acceleration g grav = GM/r 2 }156l Il57| . 

The breaking of spherical symmetry by hydrodynamic instability, manifesting itself in the growth of initially small, 
random seed perturbations to large-scale structures, turned out to be generic to the shock stagnation phase in 
collapsing stellar cores and to be decisive for the success of the neutrino-heating mechanism and for the further 
destiny of the stalled accretion shock. 


C. Neutrino-Heating Mechanism 

The development of a neutrino-heating layer is a natural consequence of the contraction of the PNS and associated 
compactification of its surrounding accretion layer during the postbounce accretion phase. The contraction of the 
PNS leads to increasing neutrinospheric temperatures and therefore growing mean energies of the radiated neutrinos 
(Fig. [7J). More energetic neutrino emission together with the decreasing postshock temperature at larger shock radii 
allows for the appearance of a gain radius: Since the temperature in the postshock layer drops roughly like r -1 
(this is well fulfilled for convectively mixed, isentropic conditions, while in ID the gradient is even steeper), the 
neutrino-cooling rate per nucleon by captures of (nondegenerate) e~ and e + on protons and neutrons drops with 
q~ oc T e oc r~ e . In contrast, the neutrino-heating rate per nucleon (which is largely dominated by v e , v e absorption 
on free n, p, respectively) scales with oc L u (e 2 )r~ 2 and thus decreases less steeply with r than q ~, enabling a 
crossing point R g to occur 0, where T^R g OC 


1. Heating Efficiency and Energetics 


With a density profile p oc r 3 between gain radius R g and shock R s (see, e.g., [l55| l and the preshock (free-fall) 
velocity vq = —\J2GM ns /R s , mass infall rate M = 4:TtR 2 |uo|po of the progenitor star, and density jump at the shock, 
/3 = pi/po ~ 10, the optical depth for v e and v e absorption in the gain layer can be estimated as 
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where Y n ss Y p ss 0.5 and (cr a bs) ~ 3.26 x 10 41 [^7^/(4 MeV)] 2 cm 2 for the average absorption cross section of a 
blackbody neutrino spectrum with temperature T v (therefore (e 2 ) ~ 21(£;bT , 1/ ) 2 ) was used. Equation Q] suggests that 
for typical accretion rates, M = 0.1-0.5Afg/s, several percent of the neutrino luminosity from the neutrinosphere can 
be absorbed in the gain layer, thus accounting for a neutrino-heating rate Q+ = r {L Vt . + L 9e ) ss 10 51 -10 52 erg/s for 
v e and v e luminosities, L v , of some 10 52 erg/s during the postbounce accretion phase (Fig. 0. 

However, in a dynamical situation as in the gain layer, where the matter is not at rest, the optical depth (which 
determines the interaction probability of a crossing neutrino) is not a perfectly appropriate measure for the heating 
efficiency. This holds in particular in the multi-dimensional case, where accretion funnels carry cold (low-entropy) 
matter from the shock towards the NS, while neutrino-heated matter expands outward in high-entropy bubbles. At 
such conditions the residence time of the matter in the gain layer accounts for the duration of its exposure to neutrino 
heating. While in the ID case the adve ction time f a d v ~ ( R s — R g )/\vi\ with v\ = vo//3 measures how long the 
accretion flow needs from R s to R g EH [158|, the dwell time in the gain region is better captured in the multi-D 
situation by the more general expression jl04 Il0(| 
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which relates the mass in the gain layer M g with the mass accretion rate M through the shock and (for conditions 
near steady state) through the gain radius. With an energy-transfer rate per nucleon by neutrino absorption of 
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q + = n !/ (e i y)(cr a bs)c, where n v Q v ) = L^/(ATrr 2 c), each nucleon absorbs an energy of qt Idweii ~ 50MeV for time 
fdweii ~ 0.1s when k B T v = 4MeV, L u = 3 x 10 52 erg/s {v £ {v e ,u e }), and r = R g ~ 100 km are assumed. This 
corresponds to a temperature of ^3 MeV and an entropy of ~20 fee per nucleon of an e -photon-dominated plasma 
at p ~ 10 9 g/cm 3 (cf. Figs. U and© • The total energy deposition rate by v e plus v e absorption in the gain layer thus 
becomes 


Qt = qt — 


9.4 x 10 51 ^ 
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k B T v \ 2 / L v \ f M g \ / R g 
4 MeV ) ^ 3 • 10 52 erg/s ) ^ 0.01 M 0 ) V100 km 
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Equation (J 6 J) corresponds to a heating efficiency of 
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and an integral energy deposition of 
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These numbers, reduced by 20-30% for neutrino-cooling losses in the gain layer, are well compatible with results of 
detailed simulations [ 8 l|, ll06( | . 


2. Hydrodynamical Explosion Models 


Neutrino-driven explosions can be found in ID simulations only for ECSNe of low-mass progenitors (EMmaiifia, 
considering an 8.8 M 0 star with an ONeMg core Because of the extremely steep density gradient at the edge 

of the ONeMg core and the corresponding rapid decrease of M, the radius of the accretion shock grows continuously 
(cf. Eq. [2]), thus creating ideal conditions for neutrino-energy transfer (Fig. [5J. The latter drives a baryonic outflow, 
which carries the energy for the explosion. With the most sophisticated available t reatment of neutrino-matter 
interactions (Table Q] Sect. IIIIBl) an explosion energy EoNeMg ~ 10 5O erg was obtained fioi UHi , which is enhanced 
at most by ~10% in 2D models due to a brief phase (~50-200ms p.b.) of convective overturn behind the rapidly 
expanding shock [3^ . 1l59( . The low explosion energy and little nickel ejection (several 10 _3 Mq; 33|) are compatible 
with estimates for the Crab SN ([111 and refs, therein). 

For more massive Fe-core progenitors nonradial hydrodynamic instabilities —convective overturn @, i m |73| in 
com bination with SASI activity [841 ]— turned out to be decisive for the success of the neutrino-heating mechanism [1041 . 
Il06| . While ID models did not explode, the Garching group found neutrino-driven, though weak, explosions for 11.2 
and 15 Mg stars in 2D simulations pi iioa. Recently, these results with the Prometheus-VERTEX program (80|, HI] 
could be_ confirmed by general relativistic 2D simulations [81] based on the newly developed CoCoNut-VERTEX 
code [l03| . which also produced explosions for solar-metallicity 27 M 0 and metal-poor (10 -4 solar metallicity) 8.1 M 0 
progenitors with Fe cores (see Fig. [ 6 ] and |l 6 l[ l. 

Neutrino-driven explosions for a variety of stars were also seen in 2D and 3D simulations of other groups with 
different multi-group treatments of neutrino transport iaang, whereas the Arizo na-J erusalem-Princeton (AJP) 
collaboration could not observe any success of the neutrino-heating mechanism paiioa. This underlines the sen¬ 
sitivity of the outcome qualitatively and quantitatively to details of the input and methods. While there are many 
differences between the modeling approaches concerning numerics and nricrophysics, the Garching 2D models in par¬ 
ticular include the full set of neutrino-matter interactions of Table [T] and a careful implementation of all fluid-motion 
dependent terms as well as GR effects in the transport. All were recognized to be significant [79], 180] but were simplified 
or ignored in the AJP calculations due to the omission of energy-bin coupling there (see also Sect. IIII Bl) . Instead of 
attributing conflicting results to transport differences, it has been repeatedly stated that the success of the Garching 
simulations is disputable because the unacceptably soft LSI 8 O-E 0 S was used (e.g., HH). This criticism, however, 
misses foundation because of the arguments given in Sect. IIII Cl and the fact that the 11.2 M 0 explosion could be 
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reproduced with the stiffer STOS-EoS (Fig. [4]), which did not lead to an explosion in the 11.2 Mq run of the AJP 
collaboration jl09j . 

Self-induced neutrino-flavor conversions in the SN core, which could swap lower-energetic u e and v e spectra with 
hotter v x and v x spectra and thus might enhance the neutrino heating behind the shock and strengthen the SN ex¬ 
plosion, have recently been shown not to have an impact during the postbounce accretion phase. Because detailed SN 
models yield electron densities that are higher than the neutrino densities (mostly even n e n v ), the matter back¬ 
ground dominates and thus suppr esses collective neutrino oscillations by dephasing the flavor evolution of neutrinos 
travelling on different trajectories [l63l - il66 |. 


3. Effects of Nonspherical Flows 


Nonradial, turbulent flows increase the residence time of matter in the gain layer HEM 1H3 and thus the mass 
M g in this region (for given M; Eq. [5]). This leads to a higher total energy deposition rate by neutrinos, Q+ (Eq. [6]), 
and to an enhanced integral energy transfer En (Eq. |5J). 

Rayleigh-Taylor fingers, for example, which develop in a convectively unstable situation (Eq. [3]) channel cool, freshly 
accreted material from immediately downstream of the shock towards the PNS and thus closer to the gain radius, 
where neutrino heating is strongest. At the same time, expanding bubbles of buoyant, high-entropy gas allow freshly 
heated matter to rise away from the gain radius instead of being accreted inward to the cooling layer. This reduces 
energy losses by the reemission of neutrinos, which can have important dynamical consequences because cooling for 
r < R g is usually much larger than net heating in the gain layer. The combination of all such effects leads to an 
increase of the temperature and pressure in the gain layer, which in turn pushes the shock farther out. A positive 
feedback cycle is the consequence, which for sufficiently strong neutrino heating enables an explosion in the multi-D 
case even when the neutrino-driven mechanism fails in ID. 

SASI activity can basically have the similar effects. It is not only associated with shock expansion and nonradial 
mass flows, thus allowing for a larger efficiency of neutrino-energy deposition, it also leads to secondary shocks that 
dissipate kinetic energy and produce e xtra heating and higher entropies, strengthening the convective activity and 
giving additional push to the shock Sim. 

After the onset of the explosion, the nonspherical situation permits simultaneous shock expansion and ongoing 
accretion. This maintains higher neutrino fluxes and stronger neutrino heating for a longer time compared to the ID 
case Il06l |. where the accretion luminosity decays as soon as shock expansion quenches the mass infall to the PNS. 

While hydrodynamical simulations clearly de mons trate that violent convective and SASI activity are crucial for 
the success of the neutrino-heating mechanism [16l| . the nature of the SASI and the exact role of hydrodynamic 
instabilities and turbulent motions for the onset of the explosion are still a matter of intense research. The SASI, 
whose amplitude grows from small pressure and/or entropy perturbations in an oscillatory way with hig hest growth 
rates for the lowest-order spherical harmonics (dipolar £ = 1 and quadrupolar £ = 2) components (84lll68i - ll70l |. seems 
to be an ideal candidate to explain not only the global asymmetries found in the SN core in simulations, but also the 
large asphericities and mixing effects that are observed in most SNe and SN remnants and that are probably linked to 
the measured high kick velocities of ma ny yo ung pulsars (cf. Sect. IVT1) . Linear growth analysis, numerical toy models 
for the linear and nonlinear regimes [76j, 17l|, and laboratory shallow -wat er experiments [172| yield evidence that the 
underlying instability mechanism is an advective-acoustic cycle fMImt rather than a purely acoustic one 


Runaway Threshold 


Burrows & Goshy [175l | conjectured the transition to the explosion to be a global instability of the postshock layer. 
Considering steady-state accretion conditions in ID, they determined a critical neutrino luminosity L V ^ C (M) as a 
function of the mass accretion rate that damps shock expansion (cf. Eq. [2j), above which they could not find accretion 
solutions and thus expected the o nset of the explosion. Their reasoning is su ppo rted by subsequent similar analyses 
of stationary accretion flows |l76l - ll78 | as well as hydrodynamical studies [TTII87L [88l Il67l . Il79j . 

The existence of a functional relation L„ jC (M) as threshold conditio n to an exp losion can be qualitatively understood 
by simple analytic considerations. Numerical simulations dMEMEM ll79L ll80( have shown that the necessary 
condition for runaway expansion is given by i a dvAheat > 1, he., the explosion can set in when the advection timescale 
of the mass flow through the gain layer, 


Gdv 




Rl /2 

V2 GM ns 1 


(9) 
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exceeds the heating timescale for neutrinos to deposit enough energy to lift matter from a gravitationally bound 
to an unbound state. In this context the total energy of the gas is the relevant quantity, i.e., the sum of internal, 
gravitational, and kinetic energies, which in a bound state is negative. Making the assumption that this energy scales 
roughly with the gravitational energy near the gain radius, which for a nucleon is e g = —GM ns m u /R g , one obtains 


Jfgl MnsRg 

qt L v {e*) ' 


( 10 ) 


Setting both timescales equal, t a dv = iheat, and using Eq. @ for R s and the fact that R g follows approximately R ns , 
i.e., R g oc R ns oc L\! 2 /{k^,T v ) 2 (which was also applied in deriving Eq. [2]), leads to 


L„, C (M) oc P~ 2,b M 2/5 M * s /5 . (11) 

This relation reproduces the functional behavior found in eh very well (the numerical factor of the scaling relation 
becomes (5-6) xlO 52 erg/s for /3 ~ 10, M = lM 0 /s, and M ns = 1.5 M@, slightly varying with the choice of other 
involved parameters). It must be pointed out, however, that the limiting luminosity for steady-state accretion solutions 
as derived in jl75j | was shown to be usually close to, but not identical with the runaway threshold at f a dv > theat- 
The latter roughly coincides with the time when the fluid behind the shock begins to develop positive total specific 
energy (see the excellent study of dz 3 ). 

Numerous stud ies f or both stationary accretion flows [M imt and time-dependent conditions in collapsing stellar 
cores 0, [§3, HE Il67t | have demonstrated that the critical neutrino luminosity for fixed mass accretion rate is signif¬ 
icantly lowered in the multi-dimensional case, typically by several 10%. The possible (or combined) reasons for this 
improvement were discussed in Sect. IIV C 31 but many aspects are still unsettled. For example, the properties and 
conseque nces of neutrino-driven turbulence (e.g., convective energy transport and pressure) are a matter of ongoing 
research |l8l| and the effects of 3D hydrodynamics have not been clarified yet. While there is hope that these might 
make the runaway easier than in 2D [83,Ha| and thus lead to earlier and mo re p owerful explosions, not all studies 
revealed a significant reduction of the threshold luminosity in 3D relative to 2D [88[. The 2D/3D comparison obviously 
depends on subtle differences of the background flow, neutrino source terms or even numerics, and requires further 
exploration. A sophisticated neutrino transport seems necessary for reliable answers. 


5. Modes of Global Instability 

The results of SN simulations and analytic studies suggest that the onset of the explosion is connected to a global 
runaway instability of the postshock accretion layer [l30l Il75j fueled by neutrino energy deposition above a certain 
threshold (see Sect. IIV G 4[) . An important question concerns the type of mode that grows fastest to trigger the 
runaway [179(. Unstable oscillatory and non-adiabatic rad ial modes were observed in many time-dependent ID 
simulations - in agreement with linear stability analysis [l82| - for neutrino luminosities intermediate between those 
that are too low to drive explosions an d those that suffice to trigger an explosion by the instability of a nonoscillatory 
radial mode 0, HE HE HE Il67l Il69l Il79l Il83j . But what happens in the multi-dimensional case? Is the runaway 
there caused by an unstable radial oscillatory or nonoscillatory mode, whose development is affected by turbulence 
altering the conditions of the background flow? Or is a nonradial nonoscillatory (possibly convective) or oscillatory 
(SASI) mode decisive? Exploration of the growth conditions has only begun, suggesting that unstable large-scale, 
nonoscillatory modes require the highest driving luminosities Il82| . but their growth may strongly de pend on the 
conditions in the SN core ML the dimensionality of the problem, and even a modest rate of rotation [184l |. While 
first 3D simulations 0, mm Hi HE HE show strongly damped or no radial oscillations, sugge sting that SASI 
modes are less strong in 3D and the explosion might be connected to unstable nonoscillatory modes [l85l Il86j | , none 
of these simulations was performed with a combination of sufficiently sophisticated neutrino transport, high enough 
numerical resolution, and a consistent inclusion of all dissipative processes (such a the decay of the NS core luminosity, 
changes of the accretion luminosity, and the shrinking of the nascent NS, all of which provide a negative feedback) 
included and combined consistently. Final answers will require well-resolved, full-scale 3D radiation-hydrodynamics 
simulations with reliable neutrino treatment. 


D. Magnetohydrodynamic Mechanisms 

MHD phenomena, in particular the magnetorotational mechanism (MRM) prop osed in EH EH , have been dis¬ 
cussed as potential drivers of SN explosions already in the 1970’s (e.g., jl89l . Il90j f and were explored by first axi- 
symmetric simulations with approximate microphysics and artificially imposed stellar core rotation and magnetic field 
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configurations in 

the neutrino physics (e.g., 19114194 and refs, therein), or more recently used neutrino transport by MGFLD (assum- 


. T hese and a flood of subsequent 2D calculations, which either ignored or radically simplified 


ing, inappropriately, the stellar medium to be at rest; EH), have confirmed that MHD processes and especially the 
MRM are able viable agents to extract energy from a highly magnetized NS and to violently expel the outer stellar 
layers. 

Because of the extremely low resistivity of SN matter, magnetic field lines are frozen in the flow. Magnetic flux 
conservation therefore leads to compressional amplification of the average strength of the nonradial field during CC, 
B oc Rcoie oc Pcore, and a corresponding growth of the energy density of the magnetic fields (oc B 2 ). Initial fields 
as expected in stellar cores at the onset of gravitational instability, i.e., several 10 9 G for the dominant toroidal 
component [36[> can thus grow by a factor >1000 but cannot gain dynamically relevant strength, for which the 
magnetic pressure has to reach a fair fraction of the gas pressure. 

Therefore secondary amplification mechanisms are crucial to further boost the magnetic energy density to values 
close to equipartition with the stellar plasma. In the MRM such an increase in energy density can be achieved 
by tapping the huge reservoir of rotational energy, E rot oc J 2 ole /(M cole R 2 ore ), that builds up at the expense of 
gravitational energy due to angular momentum (J core ) conservation during the infall. The rotational energy in a 
rapidly spinning PNS with average revolution period P ns thus becomes 


E Iot ~ 2.4 x 10 52 erg 


M n 


1.5 M 0 


Rns V 

10 km ) 


1 ms 


( 12 ) 


The amplification can either happen through the wrapping of an (initially present or convectively created) poloidal 
field, stretching it into a toroidal one, which leads to a linear increase with the number of windings. Or it can take 
place by exponential amplification with characteristic timescale of order 47r|clfl/d lnr| _1 (with D(r) = 2n/ P m t(r) 
being the angular frequency for local spin period P ro t) through the magnetorotational instability (MRI; ll95l 11961 ]), 
whose growth conditions in SN cores were studied in detail in [197j . Both processes require differential rotation, which 
naturally develops during infall. Saturation fields of order 


B 2 ~ 4t xpr 2 Vt 2 


d In fl 
d In r 


(13) 


can be expected in an MRI-unstable environment, in which dD/d lnr < 0 must hold to enable the growth of long- 
wavelength, slow-magnetosonic waves. For sufficiently large angular velocities, fields of order 10 15 -10 16 G were esti¬ 
mated [l96j . 

The ejection of matter can be driven by magnetic pressure and hoop stresses, magne tic buoyancy, or gas heating 
due to the dissipation of rotational energy through turbulent magnetic viscosity jl8dl . Il89l ll96j| . Well collimated, 
bipolar outflows or jets along the rotation axis with characteristic power 


Smhd - 10 52 f (io§g) (^) 
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10 3 rads - 


(14) 


may be generic [EH EMEU- 

Since the MRM can tap only the free energy of differential rotation in the PNS, -C E rot , angular velocities 
near the Keplerian rate of the progenitor core (P core ~ 1 s) are required for magnetic fields to grow to dynamical 
significance. SN simulations jllOl Il80j suggest that the spin period must be P core < 2-5 s, leading to NSs rotation 
periods of P ns ~ (i? n s/Rcore ) 2 Fcore under the assumption of strict angular momentum conservation. Present stellar 
evolution models that include angular momentum loss through magnetic breaking, however, yield typical core-rotation 
periods of P CO re > 100 s before collap se (c f. Sect. Ill Bit . Such slowly rotating stellar cores are consistent with obs erved 
spin rates of newly born white dwarfs 1190} and with the estimated spin periods of new-born pulsars of ~10 ms (36l . l200| . 
but they are much too slow to provide the rotational energy reservoir for powering SNe through the MRM (see Ea. fl2l) . 

A variety of mechanisms have also been proposed for magnetic field amplification in collapsing cores with no or 
slow rotation, e.g. by a convective dynamo 201 ], turbulent SASI motions in the postshock region [2021 ] or expon entia l 
steepening of Alfven waves created by fluid perturbations at Alfven points in the accretion flow of the PNS [203| . 
Moreover, Alfven waves emitted from the convective layer inside the PNS (thus extracting energy from the rich 
reservoir of gravitational binding energ y of the contracting remnant) were suggested to provide extra energy to the 
stalled SN shock by dissipative heating |204| similar to the heating of the solar corona by Alfven waves emerging from 
the surface of the Sun. Recent 2D CC simulations with neutrino transport [l!4j . however, find that these mechanisms 
are either inefficient or able to amplify the fields only locally. Large-scale fields with dynamical importance must 
reach magnetar strength (at least 10 14 -10 15 G) but in the absence of magnetorotational processes seem to require 
pre-collapse fields 100 times stronger than predicted by stellar evolution models. 
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Magnetic fields are therefore likely to play some role in all SN cores, but at the moment it seems only certain that 
they are crucial for the explosion of very rapidly spinning stars, which are probably linked to long GRBs and HNe 
(Sect. Ill Cl) . MHD mechanisms have the advantage of not being strongly coupled to the mass in the gain layer and 
the mass-accretion rate through the stalled shock, which determine the explosion energy of SNe powered by neutrino 
heating (Eq. [8j). MHD-driven explosions can therefore be considerably more energetic than neutrino-driven SNe, 
where blast-wave energies of ~(l-2)x 10 51 erg may be the upper limit (see Sect. I VI Cl) . Large global deformation and 
well collimated jets must be expected to be generic to MHD explosions of very rapidly rotating stellar cores and seem 
to be characteristic of most hyperenergetic SNe Ib/c. 

Reliable and predictive multi-dimensional simulations of such phenomena are hampered by the fact that the true 
nature of MHD phenomena can only be treated in 3D, and such models should also include reasonably realistic 
neutrino transport. Another problem arises from the extreme dependence of the dyna mical evolution on the initial 
conditions, in particular the rotation rate and profile of the stellar core (e.g., 19ll4l93| |l and the initial strength and 
geometry of the magnetic field (e.g., |l94l | and refs, therein). Moreover, a large variety of MHD instabilities, among 
them the MRI, demand high numerical resolution, which is particularly hard to achieve in 3D models, and which adds 
to the computational demands that result from long evolution times on the one hand and severe time-step constraints 
set by high Alfven speeds, v\ = B/y/Airp, on the other. The exploration of magnetorotational explosions will therefore 
remain a challenging task for the coming years. 


E. Acoustic Mechanism 


A new CCSN mechanism was envisioned in [l08l Il09l | based on results of 2D hydrodynamic simulations, which 
did not yield explosions by neutrino-energy deposition. At late times after bounce (>ls), large-amplitude dipole 
(£ = 1) gravity-mode oscillations of the PNS core were found to be excited by SASI sloshing motions of the postshock 
layer and by anisotropic accretion downdrafts. The PNS vibrations (with amplitudes of several km) were damped by 
sending strong sound waves into the surrounding medium. Running down the density gradient away from the PNS 
the waves steepened into secondary shocks. The dissipation of the latter helped to heat the postshock region. Thus 
robust explosions were obtained for a variety of progenitors. For the conversion rate of accretion power, 
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(see also eq. 1 in |l08j |b Here, R 0 is the accretion-stream radius, H 0 the wave height, pn = p/(10 n g/cm 3 ) the 
average density at the “surface” of the PNS core, and p n s,i3 = ffns/(10 13 cm/s 2 ) the average gravitational acceleration 
( g ns = GM ns /R ns ) at the PNS “surface”. The reference value of E SO und of Eq. (flfill is suggestive. It exceeded the 
neutrino-energy deposition rate (~ r (L Ve + L Pe ); Sect. IIV C ID at late times in the numerical models. The violently 
vibrating PNS thus acted as a transducer channelling accretion energy efficiently into sound. 

The fraction of the accretion power that is conver ted into core g-mode activity of the PNS could not be extracted 
reliably from the numerical calculations of [108l . Eq1, thus leaving the value of Hq uncertain. Also final numbers for 
the explosion energies could not be determined, but the 2D explosions occurred very late, implying large NS masses 
and tending to be low-energetic. Fundamental questions about the excitation efficiency of the large-amplitude, low- 
order g-modes in the PNS remain to be answered, in particular whether 3D simul ation s would yield this phenomenon 
as well. So far other groups have not been able to reproduce the results (e.g., (l0d |). maybe because their models 
either were not evolved to sufficiently late times or exploded by neutrino heating before. A serious counterargument 
to the proposed scenario was made in |2Q.6j. Employing pertubation analysis the authors concluded that non-linear 
coupling between the low-order primary modes and pairs of high-order g-modes, whose small wavelengths cannot be 
resolved in hydrodynamical simulations, damps the low-order mode amplitudes to dynamically insignificant size. The 
thermalized pulsational energy is lost by neutrino emission. 
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F. Phase-Transition Mechanism 

A first-order hadron-to-quark matter phase transition that occurs at a sufficiently low density can have dynamica l 
consequences during the postbounce accretion phase of the collapsing stellar core. This was discovered by [2071 . l208j 
using a hybrid EoS with a mixed phase that was softer than the hadronic phase and the pure quark phase. The 
latter was described by suitable choices of the parameters in the MIT bag model for strange (u,d,s) quark matter. 
Different from laboratory (heavy-ion collision) conditions, where the proton fraction Y p is close to 0.5, the mixed phase 
appears at subnuclear densities for SN matter with Y p < 0.3, and for all proton-to-baryon ratios shows a decrease 
of the transition density with higher temperatures. This is in contrast to other models for the hadro n-qua rk phase 
transition, which predict an increase of the phase-transition density with increasing temperature (e.g., [209j, |21Q|). 

The very special properties of the hybrid EoS lead to gravitational instability of the PNS when it has accreted 
enough matter and has heated up during its contraction, entering the transition to quark matter in a growing dense- 
core region. The decrease of the effective adiabatic index there below the critical value for stability triggers a second, 
supersonic implosion to the denser pure quark phase, where the EoS suddenly stiffens again. This leads to considerable 
release of gravitational binding energy and the formation of a strong, second bounce shock, which catches up with 
the stalled primary shock to cause a SN explosion even in ID models. When the second shock breaks out of the 
neutrinosphe res, e + captures by neutrons in the shock-heated matter emit a v e burst that may be detectable for a 
Galactic SN (21 1( | . 

Though this is an interesting, new scenario, the fine tuning of the QCD phase transition is problematic. In particular, 
all Eo S ve rsions that lead to explosions so far are not compatible with the 1.97±0.04Mq NS mass limit of PSR J1614- 
2230 |l24j . Changing the EoS parameters to reduce this inconsistency leads to a larger radius of the hybrid star 
and a less extreme density difference between hadronic and pure quark phases (208l |. Whether SN explosions can be 
obtained with deconfinement scenarios compatible with PSR J1614-2230 still needs to be shown. 


V. OBSERVABLE SIGNALS FROM THE SUPERNOVA CORE 

Neutrinos, gravitational waves, and heavy-element formation in the neutrino-heated ejecta are direct probes of the 
processes in the SN core. Because of the gain in sophistication of the models, in particular in the neutrino transport 
and by the inclusion of relativistic effects (recently also in 2D simulations), and because of the growing understanding 
of hydrodynamic instabilities during the postbounce accretion phase, interesting new aspects were discovered and are 
fit to change our notion of some of the accompanying phenomena and implications of CCSNe fundamentally. 


A. Neutrinos 


Neutrinos and antineutrinos of all flavors radiated from the SN core (see production processes in Table Q} carry 
information of the thermodynamic conditions (temperature, degeneracy) there, but also reflect dynamical processes 
associated with NS formation. A measurement of a neutrino signal from a future Galactic event could thus help to 
unravel the explosion mechanism. 

The shock-breakout burst of v e is a well known example of such a dynamical feature. It emerges when a flood of 
neutrinos freshly produced in shock-heated matter can suddenly escape when the bounce-shock reaches the neutrino- 
transparent regime at sufficiently low densities. Another example, though more on the exotic side, is the neutrino flash 
connected to the hadron-quark phase transition described in Sect. IIVFI Also t he p ronounced rise of the mean energy 
of heavy-lepton neutrinos, v Xl when a PNS approaches collapse to a BH [2121 . |213 | and possible —so far unclarified— 
signatures of magnetohydrodynamics can be mentioned here. 

Moreover, the large-amplitude radial oscillations [80] or low-multipole (dipolar, quadrupolar) nonradial oscillations 
of the postshock layer (due to the SASI or due to convective activity for higher-multipole asymmetries) lead to quasi- 
periodic variations of the accretion flow to t he PNS and to corresponding fluctuations of the accretion luminosity 


and mean energies of radiated neutrinos 


The effect is particularly strong for v e and i/ e , for which a significant 


part of the luminosity (amounting to a fair fraction of the progenitor-specific accretion power, Eq. 1151) originates from 
e ± captures on free nucleons in the hot accretion layer. In Fig. [3 these fluctuations are visible at t > 100 ms p.b. 
Note that the accretion luminosities depend strongly on the progenitor, and in both of the 11.2 and 15 Mq models, 
accretion continues until the end of the simulated evolution simultaneously with the accelerating expansion of the 
shock (Fig.® and with the development of positive ejecta energy. Ongoing accretion is obvious because the v e and v e 
luminosities are still c onsiderab l y hig her than those of v x , whereas after accretion has ceased all luminosities become 
nearly equal (Fig. [5] (l06l ll07l . ll60| V The “luminosities” in the left panels of Fig. [3 are energy-loss rates, E v , of 
the PNS ( v E {u e , v e , v ) j Jl u T , v T }) and not observable fluxes, the variation amplitudes are therefore damped by the 
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FIG. 7: Neutrino signals from general relativistic 2D simulations of core collapse and explosion of 11.2 Mq (upper plot) and 
15 Mq (lower plot) stars shown in Fig. [6] [l47tl . The left panels of each plot show luminosities (i.e., total neutrino-energy loss 
rates of the PNS; upper panels) and mean energies (defined by the ratio of total neutrino energy-loss rate to number-loss rate, 
E v /N v ; lower panels) with black lines for i/ e , red for v e , and blue for one kind of heavy-lepton neutrino v x . The right panels 
display the corresponding relative hemispheric differences after core bounce (the infall remains spherical). All quantities are 
measured in the lab frame at large distance. Note that the fluctuations, sudden jumps, and north-south differences at t > 300 ms 
in the upper plot are caused by violent, time-dependent, anisotropic downflows and corresponding changes of the accretion rate 
of the PNS. 


integration over all directions, and their true magnitude (percents to tens of percents) can be better read off the right 
panel, where relative hemispheric differences are displayed. 

The SASI and convective modulation of the neutrino emi ssion was not only seen in 2D simulations with RbR 
neutrino transport, but also with multi-angle treatment ( jl!2| ; cf. Sect. IIII PI for a discussion of methods). It may be 
detectable for a Galactic SN at a fiducial distance of lOkpc with IceCube or future megaton-class instruments |215| . 
First 3D simulations with approximative neutrino transport suggest that the variation amplitudes could be smaller 
than in 2D (96j j , but more and better 3D models with multi-group transport are needed for reliable information. 

Another remarkable property of th e neu trino signals in Figs. [5] and [7] is the close similarity and even crossing of 
the mean energies of v e and v x fl?)7l.l214|. Instead of the previous notion that v x are significantly more energetic 
than v e and P e , i.e., instead of a neutrino-energy hierarchy of (e„ e ) < (ep e ) < (e Ux ) with typically ( e„„) > 18-20MeV, 
state-of-the-art models now yield (e„ e ) < (ep e ) ~ (e Vx ) and (e Vx ) < 13-16MeV (Figs. [5] [3 and ref. |214j : also ref. i~7!)i . 
where rms energies are given, however). The exact value depe nds on the time and EoS: A softer EoS lead to a more 
compact and hotter PNS and thus higher mean energies |214j . 

While during the later accretion phase the order of (ep e ) and (e„ x ) can be reversed (Fig. [TJ), one obtains a mild 
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hierarchy (ep e ) < (e Vx ) during the PNS cooling phase (Fig. [5]), which, however, diminishes to insignificant differences 
at very late times (seconds after bounce), where the mean energies of all neutrino kinds become nearly identical, 
{e Ue ) ~ (ep e ) ~ (e„ x ) (Fig. 0 |l07l l216t |). The late behavior can be understood by the flat temperature profile inside 
the PNS during the late cooling stage and the close proximity of the neutrinospheric positions of all neutrinos then. 
This proximity is caused by the strong dominance of neutral-current scatterings in the effective opacity (i.e., inverse 
mean free path) for thermal coupling between neutrinos and stellar medium, K e g- = K e (n s + K e ) (n e and k s being the 
opacities for neutrino-production processes and nucleo n sca tterings, respectively) at conditions where e~ are strongly 
degenerate and neutrons start to become degenerate [216j| 4 . The close similarity of (e„ x ) and (ep B ) during the early 
postbounce and accretion phases is fostered by nucleon-nucleon-bremsstrahlung as main v x production channel |220j ], 
because it shifts the energy-sphere of u_ x to lower temperature |22l| . However, the effect is considerably enhanced 
(compare the two cases discussed in [ 10 7| ) by energy losses of v x in collisions with free nucleons N = n, p (“inelastic” 
—better “non-conservative”— nucleon recoil; Tabic HI). Such losses occur when the neutrinos diffuse out through the 
(optically) thick scattering layer between energy- and transport-spheres j22lll222l |. The small but very frequent energy 
transfers with an average value per collision of (A e v ) V N ~ £i/(6feB T — e„)/(m^rc 2 ) [223l | can force the ^-spectrum 
to become even softer than that of P e , whose production in a hot accretion layer by e + captures on neutrons is very 
efficient. 

The close similarity of the spectra of all neutrinos and the corresponding relevance of non-conservative nucleon 
recoils, which is still widely ig nored , have a bearing on the consequences of neutrino-flavor conversions, e.g., the rise 
time of the d etec table v e signal |224l| , and on neutrino-induced or flavor-oscillation dependent r-process nucleosynthesis 
in SNe [225l l226l | . Also the s teep rise of (e„ x ) before BH formation, which was found without taking non-conservative 
nucleon recoils into account [212l l213| , may be affected. 


B. Gravitational Waves 


Any nonspherical, accelerated mass motions in the SN core act as source of GWs, whose amplitude at a source 
distance D scales with the second derivative of the mass-quadruple moment, h ~ (2 G/c A )(Q/D). The GW signal 
reflects the activity phases partly visible also in the neutrino-luminosity variations (compare Figs. Q and [5j). Convective 
overturn caused by negative entropy gradients in the deceleration region of the bounce shock and in the neutrino¬ 
heating layer, the violent SASI sloshing of the whole postshock volume, and the impact of accretion dow ndrafts, 
which occur even after the onset of the explosion, induce surface g-modes in the outer layers of the PNS [214j |. whose 
frequency determines that of the GW emission [227j |. Because the buoyancy (Brunt-Vaisala) frequency connected to 
the gravity waves depends on the compactness of the PNS (see text following Eq. [3]), not only the stiffness of the EoS 
has a big influence [214| | but also GR gravity is crucial to predict the GW spectrum, which for relativistic models 
peaks around 800-1000 Hz while Newtonian simulations yield significantly lower frequencies [l47t | . 

GWs are expected to carry characteristic signatures of the explosion mechanism. While wave components associated 
with long-lasting convective and SASI activity and a broad-band signal seem typical of neutrino-driven explosions, 
the large-amplitude g-mode oscillations of the PNS core, which are essential for the acoustic mechanism, would lead 
to a dramatic rise of the GW amplitude shortly before the blast sets in [228l |. In contrast, rapid stellar core rotation 
as required for MHD explosions would produce a powerful GW burst at core bounce, possibly followed by postbounce 
emission due to triaxial instabilities [181 ]. 

The GW features and activity phases seen in 2D models are also found in 3D simulations, though the amplitudes 
may be different. Without symmetry constraints, however, the detail ed signal structure varies strongly with the 
observer direction and does not possess any template character mmmm. 


C. Heavy Elements 

Besides reprocessing shock-heated stellar layers by explosive burning, nucleosynthesis takes place in the neutrino- 
heated ejecta from the close NS vicinity. The ejecta consist of two components and have great potential for diagnostics 
of the SN mechanism. 

The first component consists of the early ejecta from the phase of shock revival. Its composition at freeze-out 
depends on the expansion timescale, which is intrinsically linked to the blast dynamics and thus to the explosion 


4 For the influence of the EoS-specific nucleon po tent ial energies in dense NS matter, which affect the /3-processes of v e and u< but were 
ignored in the models of Figs. [5] and □ see refs. |217U219| . 
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FIG. 8: Amplitudes of gravitational waves (GWs) from the general relativistic 2D simulations of core collapse and explosion of 
11.2 Mq (upper plot ) and 15 Mq (lower plot) stars shown in Fig. [6] [l47l] . The light brown lines (scaled down by factors of two 
and five in the upper and lower panel, respectively) display the growing amplitude connected with the asymmetric neutrino 
emission. The matter signal (solid black line) exhibits activity phases associated with strong, prompt postbounce convection 
(for t p b < 50 ms), increasingly violent convective and SASI mass motions in the postshock layer before the explosion sets 
in (between ~100ms and 350ms (500ms) in the 11.2 Mq (15 Mq) case), and the continued impact of asymmetric accretion 
downdrafts on the PNS after the launch of the explosion. The non-zero value of the matter signal at late times is a consequence 
of the aspherical expansion of the shocked ejecta. 


mechanism, but it also depends on the neutron-to-proton ratio set by the competition of e ± captures on nucleons and 
the inverse v e and v e captures (top two beta-processes in Table Q}. 

A good example for the relevance of these effects are recent 2D results of ONeMg-core explosions, where acceleration 
by convective buoyancy expulses early ejecta so rapidly that this material is able to retain a neutron excess sufficient 
for weak r-proces sing , in contrast to ID models where beta-reactions in the more slowly ejected plasma lift Y e close 
to 0.5 and above [331 ]. It will be interesting to explore the combination of composition and asymmetry differences of 
the early ejecta in magnetorotational explosions compared to neutrino-driven ones when eventually self-consistent, 
well resolved multi-dimensional MHD models with sophisticated neutrino transport (instead of no or highly simplified 
neutrino treatment) become available. 

The second component is the neutrino-driven wind blown off the surface of the hot PNS by neutrino-energy de¬ 
position above the neutrinosphere after the explosion has been launched. The properties of this —in the absence of 
rotation and/or strong magnetic fields— essentially spherically symmetric outflow depend on the strong gravity field 
of the NS and on the properties (luminosities and spectra) of the radiated neutrinos, which determine the strength of 
the heating Emm. Again the beta-processes of free nucleons (Table [T| set the n/p ratio of this environment. For 
sufficiently high entropy and sufficiently large neutron excess this wind may provide an interesting site for r-process 
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nucleosynthesis [233| | . 

However, besides the long-standin g pr oblem of insufficient entropy [234 l235t | . sophisticated hydrodynamic models 
find the wind to be proton-rich [l07l Il60l ]. This is a consequence of the close similarity of the spectra and luminosities 
of v e and v e d uring the PNS cooling phase, L Pe ss L Ve and (e,-> e ) ss (c Ve ) (Fig. [5] and Sect. IV AD . which enforces 
Y e > 0.5. Since (Ull 


V e 


L Ve {e Ve +2A) 


(17) 


with e v = (e^)/(e 1/ ) and A = (m„ — m p )c 2 « 1.29MeV, values of Y e < 0.5 require considerably more energetic v e 
than i/ e , satisfying £„ e — £ Ue > 4A. Recently it was shown that the nucleon potential energies in dense NS matter, 
connected to the nuclear symm etry energy, may cause sufficiently large spectral differences of v e and v e to bring the 
wind Y e slightly below 0.5 mm- But it still needs to be seen whether this reductio n allows for an r-process. 

A dominance of protons prevents r-processing but might enable a ^p-process [23fil . l237j . Active-sterile v e -v s conver¬ 
sions involving a possible sterile neutrino v s in the eV-mass range, as suggested by an anomaly of reactor Sy spectra 
and their distance and energy variation, can decrease the proton excess by removing v e and thus suppressing their 
absorption on neutrons. A recen t investigation based on an ECSN model, however, revealed only a modest effect, 
insufficient for an r-process 2381 ]. But the results depend in a complex way on the interplay between v e -v s MSW 


matter oscillations and collective vv flavor conversion, which strongly reduces the pure matter effect. More exploration 
seems necessary. 


VI. EXPLOSION PROPERTIES AND COMPACT REMNANTS 

The explosion mechanism establishes the link between progenitors stars and SNe and their remnants. It is therefore 
natural to seek for observables that might provide indirect evidence of the processes triggering the onset of the blast. 
In the following the implications and limitations of neutrino-driven explosions will be briefly discussed with respect 
to pulsar kicks, SN asymmetries, and progenitor-dependent explosion and remnant properties. 


A. Pulsar Kicks and Spins 

Young pulsars are observed to possess average space velocities around 400km/s, some of them even more than 
1000 km/s [242j [. This is too high to be understood by the breakup of binary systems in SN explosions, and natal 
kicks of the NSs are required }243l |. 

An asymmetric initiation of the explosion can naturally impart a recoil velocity to the NS due to linear momentum 
conservation. Impulsive momentum transfer by hydrodynamic forces of anisotropically expelled gas during the con¬ 
vective launch phase of the explosion, however, is hardly able to account for more than ^200 km/s even in the most 
optimistic case [2441 ] . Therefore a strong unipolar asymmetry in the progenitor star prior to collapse - in contrast to 
the higher-nrultipole asymmetries usually resulting from the stochastic and chaotic growth of hydrodynamical insta¬ 
bilities seeded by small, random perturbations (Sect. IIV C 31) has been h ypot hesized to define a preferred direction 
in which the SN blast could break out with the highest velocities [2451 I246[ . However, stellar evolution models, 
self-consistently evolved in 3D through all advanced burning stages up to gravitational instability, which could lend 
convincing support to the existence of such global, low-multipole pre-collapse asymmetries, do not exist yet. 

Anisotropic neutrino emission is another potential NS kick mechanism by which the NS could achieve a recoil 
velocity of u ns ~ £„ • 33000km/s {E v /?> x 10 53 erg)(M ns /1.5M 0 ) _1 . An asymmetry of 1% of the total neutrino 
energy loss would thus kick the NS to more than 300 km/s. The asymmetric neutrino emission ass ociat ed both with 
postbounce accretion |94], [95| and with the convective activity during the PNS cooling evolution |l57| . however, is 
highly time variable and nonstationary in space and time; therefore, it can hardly account for recoil velocities of more 
than 10 km/s. Exotic mechanisms, invoking ultrastrong NS magnetic dipole fields and non-standard, still uncertain 
neutrino properties have therefore been suggested as speculative solution (e.g., [247| and refs, therein). 

Probably the most plausible origin of the NS velocities was proposed on t he b asis of 2D explosion simulations by 
Scheck et al. (94j , l248j] . whose results received recent confirmation by 2D [2491 l250t | and 3D models [ffil l239j | . Scheck et 
al. showed that the asymmetric expulsion of gas does not only exert “contact forces” during a few 100 ms, in which the 
explosion is lauched and ejecta and PNS interact hydrodynamically, but that the asymmetric ejecta exert a long-time, 
anisotropic gravitational pull, which can accelerate the PNS over seconds to velocities of many 100 km/s (Fig. [9]). 
For particularly large aspliericity of the ejecta a NS velocity of v ns > 1000km/s was obtained j94j. A hemispheric 
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FIG. 9: NS kicks and anisotropic Ni ejection for asymmetric SN explosions in 3D simulations [ 95 I . l239l |. Upper panels: Entropy 
isosurfaces of SN shock and convective bubbles {left) and ray-casting image of the density ( middle) at t = 1.3 s after bounce. 
The deformed boundary is the shock, the viewing direction is normal to the plane of NS kick and spin vectors (white and 
black arrows), which define the plane of the entropy distribution (right). The NS (black cross) is clearly displaced from the 
geometrical center of the expanding shock towards the side of weaker explosion. It is accelerated mainly by the asymmetric 
gravitational attraction of less rapidly expanding, dense ejecta clumps (intense reddish and bluish in middle image). Left lower 
panel: Recoil velocity of the NS vs. time for four 3D explosion simulations of different stars. The acceleration continues even 
later than 3 s and kicks of >600km/s are reached. Right lower panel: Anisotropic production of radioactive j6 Ni by explosive 
nuclear burning behind the expanding shock. For large NS kicks nickel is ejected preferentially in the direction where the shock 
is stronger, i.e., opposite to the NS motion (red arrow). 


asymmetry of the mass distribution of only Am = ±10 -3 Mq in a shell expanding away from the NS from an initial 
radius r\ = 100 km with v s = 3000 kms -1 can tug the NS to a velocity of u ns ss 2GAm/(riV s ) ~ 900kms _1 95. 231) 1. 

Gravitational forces of anisotropically ejected gas can thus mediate an efficient, long-lasting acceleration of the NS, 
transferring momentum from the anisotropically ejected matter to the compact remnant. Since the NS is gravitation¬ 
ally pulled by the slower, usually denser ejecta associated with a weaker explosion shock, Wongwathanarat et al. [239| 
expect the bulk of the iron-group nuclei and of other elements heavier than 28 Si, which are explosively produced in 
the shock-heated ejecta, to be expulsed preferentially in the direction opposite to the NS motion. They predict a very 
strong asymmetry of the nickel ejection in the case of large NS kicks (see Fig. [9]), which could be an observationally 
accessible, characteristic feature of the hydrodynamical-gravitational kick mechanism. 

Asymmetrical convective downdrafts and rising bubbles as well as violent, low-multipole SASI sloshing modes, 
which have spiral components in 3D, can establish angular momentum separation between PNS and ejecta and thus 
may cause considerable PNS rotation even if the stellar core did not rotate before collapse ([82, 83|; because of the 
use of an inner boundary condition, however, these results were questioned by 251j ). Naturally, any anisotropic mass 
infall that hits the accretor not exactly head-on can exert a torque and spin up the PNS. A mass Am = 1CP 3 Mq 
that has an impact velocity Hi mp ~ y/2 GM ns /R ns ~ 10 lo cm/s and an impact parameter d = (R ns ~ 30 km when 
colliding with the NS transfers an angular momentum of AJ ns = Amvi mp d ~ 6 x 10 46 gcnr 2 /s, corresponding to a 
NS spin period of T ns = 27r/ ns /AJ ns ~ 0.2s for a typical value of the NS moment of inertia of I ns ~ 2 x 10 45 gcm 2 . 
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FIG. 10: Left panel: Asymmetric shock front (outer bluish, nearly transparent surface) and mushroom-like, high-entropy 
bubbles of neut rino-heated plasma around the central NS (dark grey surface near the middle) at 0.5 s after bounce in a 3D SN 
simulation [240 ]. The shock front has a diameter of ~4000km. An octant is cut out to show the entropy distribution (color- 
coded between ~10 and 21 fcn from blue to yellow to red) in the expanding Rayleigh-Taylor mushrooms surrounded by cooler 
accretion downdrafts. All visible structures have grown from tiny, random seed perturbations by hydrodynamic instabilities. 
Right panel: Asy mmet ric ejection of different chemical elements during the explosion of the left image, but ~9000 s later into 
the SN evolution |24Q|1 . The side length of the displayed volume is about 7.5 x 10 7 km. The largest bubbles on the left image 
have seeded the growth of the most prominent Rayleigh-Taylor fingers in the right picture, which expand with up to 4500 km/s. 
They are surrounded by the helium and hydrogen of the outer stellar shells (not visible). Together with the smaller features 
they thus carry heavier chemical elements from deep stellar layers far into more slowly expanding, lighter SN material. Blue 
filaments contain dominantly nickel, red fingers mostly oxygen, and green is associated with carbon. A mix of nickel and oxygen 
appears in pink. The whitish glow results from a contamination with other colors as a consequence of the volume rendering for 
the visualization. 


Indeed, 3D explosion simulations yield T ns in the range of hundreds of milliseconds to seconds (95ll239j . Nevertheless, 
angular momentum transferred to the PNS by hydrodynamical flows during the development of the explosion and in 
the post-explosion accretion phase is unlikely to be sufficient to account for th e est imated NS birth spin periods of 
order ~10ms, which seems to require rotation of the collapsing stellar core (3(1 l252j . Explaining a possible spin-kick 
correlation of observed NSs remains a challenge for any discussed kick mechanism connected to explosion asymmetries 
of progenitor stars with or without rotation. 


B. Supernova Asymmetries 


The large asymmetries imprinted on the ejecta by the violent, nonradial mass motions in the SN core, which precede 
and accompany the neutrino-driven revival of the blast wa ve, se ed the growth of secondary Rayleigh-Taylor instability 
in the shock-accelerated outer shells of the exploding star [2531 ]. Since the developing Rayleigh-Taylor mushrooms are 
denser than the surrounding gas, they are less decelerated than their environment and can penetrate the composition 
interfaces of the progenitor, retaining high velocities as the SN ejecta expand. Thus they carry freshly synthesized 
radioactive nickel and other heavy elements from the vicinity of the nascent NS into the outer stellar layers. Significant 
amounts of the initially innermost ejecta can be mixed deep into the helium shell and even the hydrogen layer of the 
disrupted star [254| . destroying the well-stratified onion-shell structure of the progenitor. 

In 3D simulations of a SN 1987A progenitor model, large nickel-dominated clumps (containing up to sever al 10 ~ 3 Mq 
of 56 Ni) were found to speed through the stellar hydrogen envelope with up to 4500 km/s (Fig. [TO] [24(1 ). This can 
explain mixing phenomena and asymmetries observed in SN 1987A, e.g. the detection of X-rays and 7 -rays from the 
radioactive nickel decay much earlier than predicted by ID explosion models [69j. The outward mixing of radioactive 
nickel and inward displacement of hydrogen can well account for the shape and width of the lightcurve maximum of 
SN 1987A (V. Utrobin, private communication). 

However, it is still unclear whether explosion asymmetries associated with the development of hydrodynamic insta¬ 
bilities in the SN core and the subsequent growth of mixing instabilities in the stellar envelope are able to explain the 
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FIG. 11: Explosion and remnant properties predicted by parametrized ID neutrino-driven SN simulations I24ll1 of a large 
set of progenitor stars from [sjjj]. Explosion energy (top left), time of onset (top right), baryonic remnant mass (middle left), 
neutrino-energy release by the compact remnant (middle right), and ejected Ni mass (bottom left) are shown as functions of 
stellar birth (ZAMS) mass. The lower right plot gives the compact remnant mass vs. the enclosed mass at the base of the 
O-burning shell of the progenitor. Neutrino cooling of the dense NS core was prescribed such that the properties of SN 1987A 
were roughly reproduced for ~20 Mq stars (green histogram bar). Accretion neutrino luminosity was self-consistently computed 
by approximate neutrino transport. The ticks in some panels mark masses where computed models did not explode. Bars 
of remnant masses reaching to the upper panel edge (3 Mq) and arrows in the bottom right panel signal formation of a BH 
containing the whole mass of the star at collapse. The only exception is the 37 Mq progenitor, where the explosion ejects 
~3.2Af@ while fallback creates a BH with 6.5 Mq. Blue segments indicate fallback masses and orange segments Ni-mass 
uncertainties because of unclear Ni abundance in the iz-heated ejecta. 


prolate shape of the SN 1987A ejecta cloud. It will also have to be seen whether they can account for the extremely 
fast “jet” structures observed ahead of the explosion shock in the Cassiopeia A SN remnant. Moreover, the global 
asphericity of most SNe Ib/c might require larger nonradial deformation than the asymmetric structures that can 
stochastically grow from initially small random perturbations. 


C. Neutron Stars and Black Holes 

Remnant masses and explosion properties (energy, ejected 56 Ni mass) and their systematics with the progenitor 
mass also carry informat i on on the explosion mechanism. The observational basis of determined or constrained NS 
and BH masses mmm, SN-progenitor connections [261127^. and estimated explosion parameters (e.g., p3.[25j: 
Fig. [3]) is rapidly growing. 

From the observed mass distribution of compact remnants and its possible gap between ~2 and 5 Mq at the 
boundary between NSs and BHs, it was inferred that the SN engine must launch the (neu trino-powered) explosion 
within 100-200 ms after bounce in order not to overproduce remnants in the gap | |257l . |258| . This was considered as 
argument that the mechanism is supported by Rayleigh-Taylor (convective) rather than SASI instability. However, in 
the SN core both of these nonradial instabilities occur simultaneously (7|J Il6l| and cannot be separated just on the 
basis of a timescale argument. Moreover, the population evolution models of 12571 l258j used very simple theoretical 
considerations to determine the explosion energy for early and late explosions and to estimate the fallback mass of 
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matter that initially moves outward but ultimately fails to escape because of insufficient blast-wave energy. The 
analytic theory ignores, for example, dynamical effects and the nonnegligible additional power carried by the early 
neutrino-driven wind (cf. [§5j])- Other approaches to predict mass distributions of NSs and/or BHs were based either 
on piston-driven explosions with predefined mass cut and explosion energy (e.g., (259[), or on a single-parameter 
criter ion t o distinguish progenitors that are likely to explode or not (260l |. 

In |24l| an alternative approach was adopted. Hydrodynamical simulations in ID were performed for a large set 
(roughly 100) of solar-metallicity progenitors of [22j using an analytic, time-dependent two-zone model of the cooling, 
contracting PNS, whose free parameters were calibrated such that the explosion energy and 56 Ni mass of SN 1987A 
were reproduced for stars with ZAMS mass around 20 M 0 . The effects of accretion luminosity were taken into account 
by simplified neutrino transport [H] . With this prescription all stellar collapses and possible explosions were simulated 
for at least 15 s beyond core bounce and were followed after the PNS cooling for hours to days later until the fallback 
mass was determined. 

Results of the calculations are shown in Fig. [Tl] and reveal a number of interesting insights, which, of course, depend 
on the considered progenitor set: 


• Because the stellar structure varies nonmonotonically, the SN properties depend on the progenitor mass in a 
complex way. Large differences of the explosion characteristics are possible for small mass differences. 

• Failed explosions with BH formation seem possible for progenitors below 20 M 0 , and successful SNe with NS 
formation are found also between 20 and 40 Mq. 

• Neutrino-driven explosions with energies in excess of 2 x 10 51 erg and 56 Ni production of significantly more than 
~0.1 Mq seem unlikely. 

• The time of the onset of the SN blast (measured by the moment the shock passes 500 km) varies between ~0.1 s 
and 1.1s, so it includes “early” and “late” cases. Later explosions tend to be less energetic because less mass is 
available for being heated by neutrinos. 

• The NS baryonic masses are in the range of ~1.2-2M 0 . The smallest BH, formed by fallback, contains 6.5 Mq, 
all other BHs originate from failed explosions and contain all the mass of the progenitor at collapse (>8.5M 0 ). 
The possible gap of the observed remnant distribution is clearly reproduced. 


Fallback is larger for the lower-mass progenitors where an extended hydrogen envelope leads to a stronger reverse 
shock. The result of little fallback in solar -met allicity progenitors is compatible with conclusions drawn from 
an analysis of observed double NS systems 


261 ]. 


• Although the remnant mass is an almost monotonic function of the enclosed mass at the base of the oxygen¬ 
burning shell, the latter is no reliable indicator for the fate of the star because some models with relatively small 
Si-cores do not explode. 


• Neutrino-driven explosions are fostered by big “jumps” in the stellar density and entropy profiles (cf. Fig. O, 
reducing the mass-infall rate (and ram pressure) and allowing the shock to expand (cf. Eq. [2]) . 


Certainly these results are based on ID simulations and many approximations were made. Therefore they can be only 
a very first step, but nevertheless are enlightening concerning the implications of neutrino-powered explosions. They 
challenge a number of paradigms for the progenitor-explosion and progenitor-remnant connections. In particular 
the limited blast-wave energy and nickel production support arguments in favor of another explosion mechanism 
for HNe. These events are likely to be triggered by magnetorotational processes. More research, observationally and 
theoretically, will have to clarify whether there is a continuous transition between both, associated with a varied degree 
of progenitor rotation and leading to a continuous spectrum of explosion energies that reach from the neutrino-powered 
regime of E exp <2x 10 51 erg to the hyperenergetic regime of E ex p > 10 52 erg as suggested by some phemomenological 
studies (see Fig. 0. 


VII. SUMMARY, CONCLUSIONS, OUTLOOK 

Supernova theory has made remarkable progress over the past decade, promoted by common interests of the astro-, 
particle (neutrino), nuclear, and gravitational physics communities and by an increasing number of active (young) 
researchers in the field. A deeper understanding of the physical mechanisms that initiate and fuel SN and HN 
explosions of massive stars is of crucial importance not only for establishing the progenitor-remnant connection but 
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also for predicting the properties of stellar explosions, their nucleosynthetic output, and the characteristics of their 
gravitational-wave and neutrino signals. 

The most sophisticated present simulations demonstate that neutrino-energy deposition can power ECSNe (even in 
spherical models) of ~9 Mq stars with ONeMg-cores near the lower mass limit for SN progenitors (Fig. [5]). Overall, 
the features of such explosions, e.g., low energy and little nickel production, seem to be compatible with observational 
candidates like the Crab SN and some faint transients. Multi-dimensional simulations suggest these explosions to be 
potential sources of light r-process nuclei up to silver and palladium (Sect. IV Cl) . Several groups have also reported 
successful neutrino-driven explosions (with multi-group neutrino transport) for Fe-core progenitors above 10 Mq 
(Sect. IIVC 21 Figs. 01 0. Ultimate confirmation of the viability of this mechanism for a wider range of progenitor 
masses therefore seems to be in reach. 

The onset of the explosion can be understood as a global runaway instability of the accretion layer, whose initiation 
depends on the power of neutrino-energy deposition. While the exact mode of the runaway is still a matter of 
exploration and debate (e.g., low-multipole SASI or higher-multipole convective, oscillatory or nonoscillatory?), its 
threshold in terms of the driving neutrino luminosity is lowered by nonradial fluid motions in the neutrino-heating 
layer. Such flows play a supportive role because they stretch the residence time of matter in the gain region and 
thus decrease the heating timescale and increase the efficiency of neutrino-energy deposition, leading to successful 
explosions even when sophisticated spherical models fail (Sect. IIV C 31) . The efficiency of neutrino-energy transfer, 
the growth conditions and growth rates of different hydrodynamic instabilities, and the critical luminosity threshold 
for an explosion may not only depend on the dimension and thus will ultimately require simulations in 3D, but have 
been shown to depend also on putative “details” of the physics ingredients like approximations for the energy and 
velocity dependence of the neutrino transport, the neutrino-interaction rates, general relativity, and the contraction 
of the nascent NS in response to the nuclear EoS (Sect. [TTT] Fig. 0]). Moreover, the outcome of the complex neutrino- 
hydrodynamical simulations can be sensitive to the numerical resolution, which naturally is subject to limitations in 
full-scale, multi-dimensional SN-core models. 

While detailed modeling of the processes in collapsing stars now pushes forward from the second to the third 
dimension, facing considerable computational challenges and demands mainly for the neutrino transport, a growing 
host of studies begins to explore the observational consequences of neutrino-driven explosions. In view of existing and 
upcoming big detection facilities, in particular neutrino and gravitational-wave signals (Figs. [5] [TJ [8]) are of relevance for 
SN-core diagnostics targeting a future Galactic SN. The former even have the potential to yield valuable information 
on pa rticle propertie s of the neutrinos provided the characteristics of the SN emission are sufficiently well understood 
(e.g., [224 l262H264l p. Sophisticated neutrino transport and interac tion treatments have revealed interesting signal 
features like an amazing robustness of the neutronization v e burst [262| . characteristic differences of the rise time 
of the v e and v x emission after bounce [224 |. luminosity variations associated with nonsteady flows in the accretion 
layer [l!2l l214j , and a close similarity of the lumi nosit i es an d spectra of neutrinos and antineutrinos of all flavors 
during the PNS cooling phase (Sect. IV Al Fig. 0 |l07i l216l |f with important consequences for SN nucleosynthesis 
(Sect. [VS]). 

While a Galactic SN in the near future is a realistic possibility, it will be a unique event and might not provide 
evidence of wider validity. Photometric and spectroscopic diagnostics of extragalactic SNe and of gaseous, young 
SN remnants, which reveal information on explosion energies, 56 Ni production, ejecta masses, asymmetries, and 
composition, as well as progenitor constraints (cf. Fig. 0 are therefore extremely valuable, and more is desirable. 
First-principle explosion models begin to become mature enough to be linked to such observations, a possibility that 
defines a fruitful territory for future reseach. Neutrino-driven explosion models also begin to allow for predictions of 
compact remnant (NSs and BHs) masses, kicks, and spins. 

Nonradial hydrodynamic instabilities in the collapsing stellar core, which can grow from small, random initial per¬ 
turbations before neutrino heating revives the stalled shock, lead to low-multipole asymmetries that trigger anisotropic 
and inhomogeneous expulsion of matter. Hydrodynamic instabilities in the SN core therefore do not only yield a nat¬ 
ural explanation of the origin of pulsar kicks up to more than 1000 km/s ( Sect. fVTAl Fig.0; they also seed large-scale 
mixing processes in the exploding star, accounting for the penetration of high-velocity clumps of inner-core material 
into the hydrogen and helium ejecta of well observed SN explosions fSect. IVlBl Fig. flOl). 

First results of a systematic exploration of the progenitor-supernova connection based on the neutrino-heating 
mechanism show strong sensitivity of the explosion properties on the stellar structure and, for the employed set of 
stellar models [22j, large variations even within narrow progenitor-mass intervals ISect. fVTCl Fig. ITT1). The explosion 
models can reproduce fundamental properties of the empirical remnant-mass distribution but reveal that neutrino- 
driven explosions are unlikely to explain SN energies above ~ 2 x 10 51 erg and nickel masses significantly higher than 
0.1 Mq. This underlines the need for an alternative engine that powers stellar blast waves with energies from several 
10 51 erg up to more than 10 52 erg. Such hyperenergetic events, which typically also exhibit unusually large nickel 
ejection (Fig. 0 and deformation, are most probably energized by magnetorotational effects. 

Many questions remain to be answered in this context and require more observations and theoretical work. What 
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discriminates progenitors of “normal” SNe from those of HNe? Is rapid rotation of the progenitors the crucial 
parameter? Is it connected to binary evolution? Is there a continuous spectrum of stellar explosions connecting the 
SN and HN regimes? Is a mixed mechanism, neutrino-heating in combination with magnetorotational energy transfer, 
at work in such events? 

On the theory side the mission of clarifying the SN engines is severely handicapped by the unavailability of multi¬ 
dimensional stellar evolution models with the quality to reduce the major uncertainties of the stellar structure, 
rotation, and magnetic fields at the onset of core collapse. It is clear that reliable theoretical predictions of the 
progenitor-remnant connection and of explosion properties —energies, nucleosynthetic yields, asymmetries, remnant 
masses, and neutrino and GW signals— heavily depend on a firm knowledge of the stellar conditions at the time the 
gravitational instability is reached. 
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